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Sonic  Nozzles  for  Mass  Flow  Measurement  and 
Reference  Nozzles  for  Thrust  Verification 

(AGARD  AR-321) 


Executive  Summary 


Accurate  measurement  of  massflow  and  thrust  is  essential  to  the  success  of  windtunnel  tests  supporting 
engine-airframe  aerodynamic  integration  studies.  Among  other  benefits,  optimising  engine  integration 
results  into  reductions  of  the  cruise  drag  of  an  aircraft  by  at  least  several  percent.  Hence,  mastering  this 
technique  at  an  early  stage  of  a  project  allows,  within  a  given  set  of  specification  (e.g.  the  range),  to 
save  on  the  mass  of  the  aircraft  and  consequently  on  its  cost  -  more  generally  it  contributes  to  reducing 
technical/economical  risk. 

Refinement  of  experimental  techniques  contributing  to  cruise  drag  prediction,  and  possibilities  offered 
by  small  scale  engine  simulators  are  today  somewhat  limited  by  the  aceuracy  of  massflow  and  thrust 
measurements  on  reference  nozzles,  which  have  to  be  used  at  various  stages  of  the  experiments.  Indeed, 
determining  thrust-drag  balance  with  an  accuracy  better  than  one  percent  often  requires  subtracting 
large  quantities  which  need  to  be  known  individually  within  a  few  thousandths.  This  is  today  still 
difficult  to  do,  even  in  the  simple  case  of  reference  nozzles.  Moreover,  within  cooperative  programmes, 
systematic  interfacility  bias  resulting  from  slight  differences  in  test  methodology  can  raise  complex 
issues  for  the  partners. 

For  these  reasons,  and  although  the  topic  is  far  from  being  new,  the  Fluid  Dynamics  Panel  decided  in 
1993  to  create  a  Working  Group  (WG  19)  to  report  on  the  state-of-the-art  and  make  practical 
recommendations. 

Progress  has  been  made  in  the  physical  understanding  of  some  flow  phenomena  and  a  consensus  has 
been  reached  on  how  to  proceed,  while  keeping  in  mind  usual  wind  tunnel  constraints  and  cost- 
effectiveness.  As  regards  practical  results,  it  can  be  said  that  measuring  gaseous  mass  flows  within 
±  0.1%  or  better  is  still  very  difficult.  For  most  tests,  however,  with  reasonable  care,  bias  and  random 
(repeatability)  errors  ean  be  kept  within  ±0.1%  each.  For  thrust  measurements,  these  values  must 
typically  be  doubled. 

X.  Bouis, 
Chair,  WG  19 
January  1997 


Mesures  de  debit  par  cols  soniques  et  verification  de 
la  poussee  par  tuyeres  etalon 

(AGARD  AR-321) 

Synthese 

La  mesure  precise  des  debits  et  des  poussees  contribue  de  maniere  essentielle  au  succes  des  etudes  en 
soufflerie  de  1’ integration  aerodynamique  cellule-propulsion.  Optimiser  cette  integration  permet  par 
exemple  de  reduire  de  plusieurs  pour-cent  la  trainee  aerodynamique  d’un  aeronef  en  croisiere.  La 
maitrise  de  telles  techniques  au  cours  des  premieres  etapes  d’un  projet  permet,  a  specifications  donnees 
(ex  :  le  rayon  d’ action)  de  reduire  la  masse  de  1’ avion  et  par  consequent  son  cout ;  plus  generalement, 
cette  maitrise  contribue  a  la  minimisation  du  risque  technico-economique. 

L  affinement  des  techniques  experimentales  concourant  a  la  prevision  de  la  trainee  en  croisiere,  et  les 
possibilites  offertes  par  les  simulateurs  de  moteurs  a  echelle  reduite  butent  aujourd’hui  sur  la  precision 
des  mesures  de  debit  et  de  poussee  sur  les  tuyeres  etalon  qui  interviennent  a  divers  stades  des 
experiences.  En  effet,  la  determination  d’un  bilan  poussee-tramee  avec  une  precision  superieure  au 
pour-cent  passe  en  pratique  par  la  soustraction  de  grandeurs  importantes  dont  chacune  doit  etre  connue 
a  quelques  milliemes  pres.  Ceci  est  encore  difficile  aujourd’hui,  meme  pour  le  cas  de  simples  tuyeres 
etalon.  De  plus,  dans  les  programmes  menes  en  cooperation,  les  ecarts  systematiques  entre  installations 
d’essais  resultant  de  petites  differences  de  methodologie  ne  vont  pas  sans  poser  des  problemes  aux 
partenaires. 

Pour  ces  raisons,  et  bien  que  le  sujet  ne  date  pas  d’hier,  la  Commission  de  Dynamique  des  Fluides  a 
decide  en  1993  de  demander  a  un  Groupe  de  Travail  (WG  19)  de  faire  le  point  sur  I’etat  de  Part  dans  ce 
domaine  et  d’etablir  des  recommandations  pratiques. 

Des  progres  ont  ete  accomplis  sur  la  comprehension  physique  des  phenomenes  et  un  consensus  a  ete 
obtenu  sur  les  procedures  a  suivre  compte-tenu  des  contraintes  habituelles  des  souffleries  et  du  souci 
d’un  bon  rapport  cout/efficacite.  En  pratique,  on  pent  dire  que  mesurer  des  debits  gazeux  a  mieux  que 
±0,1%  demeure  tres  difficile.  Pour  la  plupart  des  essais  cependant,  des  precautions  raisonnables 
permettent  de  contenir  I’erreur  systematique  et  I’erreur  aleatoire  (repetabilite)  chacune  dans  la  limite  de 
±0,1%.  En  ce  qui  conceme  la  mesure  de  la  poussee,  il  faut  approximativement  doubler  les  valeurs  ci- 
dessus. 


X.  Bouis, 
President  du  WG  19 
Janvier  1997 
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1.  INTRODUCTION 


1.1  Background,  objectives,  scope  of  work 


For  drag  measurement  of  aircraft  models  in  wind  tunnels,  a  precision  of  1  "drag  count", 

ACd  =  ±  1.10-^,  which  represents  about  0.3  %  of  the  total  drag  at  cruise  conditions,  is  today 
required,  and  can  actually  be  obtained,  at  least  in  repeatability  [1],  [2]. 

AGARD  Advisory  Report  184  [3]  presents  the  wind  tunnel  flow  quality  and  data  accuracy 
requirements  which  are  related  to  that  precision. 

On  models  equipped  with  engine  simulators,  and  used  to  study  engine  installation  drag,  if  the  global 
forces  on  the  model  are  measured  with  an  accuracy  of  0.3%,  the  homogeneous  accuracy  to  be  insured 
on  the  evaluation  of  the  engine  simulator  thrust  is  also  0.3%. 

The  thrust  of  an  engine  simulator  is  derived  from  a  preliminary  measurement  of  the  exhaust  flow 
thrust  on  a  static  test  bench. 

At  high  subsonic  cruise  Mach  number,  a  ratio  of  1.4  between  the  exhaust  speed  of  a  fan  flow  to  the 
flight  speed  is  representative,  so  that  an  error  of  0.3  %  on  the  net  thrust  is  equivalent  to  an  error  of 
0.3%  X  (1.4  -  1)/1.4  =  0,086  %  on  the  exhaust  flow  gross  thrust,  which  is  only  measured  at  static. 
In  fact,  taking  also  into  account  the  thrust  of  the  core  flow  of  a  turbofan  exhaust  configuration,  it  is 
the  near  value  of  0.1%  which  can  be  retained  as  an  objective  for  thrust  measurement  on  a  static  test 
bench. 

Consistently,  the  same  precision  of  0.1%  should  also  be  aimed  at  on  any  static  test  bench  for  turbofan 
exhaust  models. 

Since  thrust  is  generally  expressed  as  a  ratio  of  thrust  to  mass  flow,  the  precision  of  0.1  %  equally 
applies  to  mass  flow  measurement. 

Mass  flow  is  usually  measured  by  the  way  of  a  sonic  throat,  and  thrust  test  benches  are  often  checked 
by  testing  a  reference  nozzle. 

The  accuracies  of  these  measurements  are  assessed  by  theoretical  calculations  and/or  by  calibration 
process. 

As  different  sonic  throat  profiles,  different  flow  calculation  methods,  and  different  calibration 
techniques  are  used  among  the  various  NATO  countries,  a  working  group  was  initiated  by  the  Fluid 
Dynamic  Panel  of  AGARD  to  present  the  state  of  the  art,  and  to  compare  different  possibilities  and 
attainable  precisions. 

It  is  the  subject  of  this  report,  which  includes  the  next  chapters  : 

2  -  Basic  calibration 

2. 1  -  Mass-time  calibration  of  massflowmeters 

2.2  -  Calibration  of  massflowmeters  by  using  transfer- standard  nozzles 

2.3  -  Nozzle  thrust  measuring  benches 

2.4  -  General  precautions  for  all  tests 

3  -  Flow  analysis 

3.1  -  Cylindrical  throat  nozzles  :  the  ASME  Long  Radius  Nozzle 

3.2  -  Toroidal  throat  nozzles 

3.3  -  Real  gas  effects 

3.4  -  Viscous  effects  on  thrust 

3.5  -  Practical  formulce  to  calculate  mass  flow  rate  and  thrust 
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4  -  Examples  of  experimental  results 

4.1-  Boeing's  experience  on  ASME  and  "cubic"  nozzles 

4.2-  3"  "cubic"  nozzle  test  results  of  different  calibration  tanks 

4.3-  Discussion 

5  -  Conclusions 


12  General  remarks 


1-  Reference  should  be  made,  before  starting  the  presentation,  to  the  existing  international  standards 
on  mass  flow  measurements  (ISO-4006,  ISO-5167,  ISO-5168,  ISO-9300),  and  the  several  excellent 
review  papers  ([4],  [23],  [25]),  some  of  them  being  nearly  30  years  old. 

2-  It  should  be  stressed  that  the  goal  of  mass  flow  and  thrust  measurements  with  an  accuracy  of 
±0.1%  is  a  serious  chdlenge  and  that  careful  application  of  ISO  standards  is  necessary  but  far  from 
being  sufficient  to  get  it.  ISO  generally  refers  to  errors  in  the  range  of  0.5%  to  2%. 

3-  A  repeatability  of  the  order  of  0.1%  to  0.2%  is  however  currently  achieved  in  some  wind  tunnels 
which  can  offer  with  considerable  precautions  a  level  of  drag  measurement  repeatability  close  to  1 
drag  count.  This  indicates  that  if  systematic  errors  are  properly  eliminated,  the  above  serious  challenge 
is  not  a  dream. 

4-  Above-mentioned  articles  should  already  be  known  by  the  reader  of  this  AGARD  document.  They 
are  indeed  the  starting  point  for  the  following  discussions  and  most  of  their  analysis, 
recommendations  and  conclusions  are  still  valid  today. 

5-  The  case  of  cryogenic  wind  tunnels  is  not  adressed  in  the  report  although  engine/jet  simulation  is 
planned  in  some  of  these  facilities,  e.g.  ETW  and  KKK.  It  poses  some  additional  problems  which 
could  not  be  fully  identified  within  the  scope  of  this  working  group. 
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2.  BASIC  CALIBRATION 


Basic  calibration  for  a  massflowmeter  means  that  direct  measurements  of  mass  and  time  are  used  to 
determine  the  mass  flow  rate.  A  number  of  facilities  over  the  world  are  able  to  provide  such 
calibrations,  mostly  for  the  natural  gas  industries,  but  very  few  with  air  as  fluid  and  giving  the 
required  acci^acy.  The  working  group  could  only  find  two  such  facilities,  one  at  the  Colorado 
p^neenng  Expenment  Station  Inc.,  CEESI,  in  the  USA  and  the  other  at  the  National  Engineering 
Laboratory,  NEL  m  the  UK.  Both  are  briefly  described  hereafter.  The  mass  flow  ranges  of  these 
facilities  are  limited  and  it  is  impossible  to  calibrate  each  massflowmeter  in  such  a  facility.  So  transfer 
standard  nozzles,  once  a  basic  calibration  has  been  performed,  are  used  as  references.  Although 
simple  combinations  of  nozzles  in  series  and/or  in  pardlel  do  give  proper  answers  if  repeatability  only 
is  looked  at,  the  target  of  0.1%  in  absolute  accuracy  can  only  be  reached  if  a  sophisticated 
methodology  is  applied,  like  that  of  Boeing  described  in  2.2. 


2.1  Mass-time  calibration  of  massflowmeters 


2.1.1  CEESI  gravimetric  calibration  facility 

This  facility  currently  seems  to  have  the  most  accurate  system  to  calibrate  venturis  or  nozzles  with  air 
It  IS  presented  in  line  diagram  form  in  figure  6. 

The  principle  is  to  discharge  air  from  a  tank  over  a  measured  time  interval  at  constant  conditions 
through  the  ventun  or  nozzle  under  test.  The  tank  is  accurately  weighed  after  filling,  and  again  after 
air  has  hem  discharged  to  determine  the  mass  of  air  used  during  the  time  interval.  Details  are 
presented  in  ref.  [8]. 

Mass  flow  is  limited  to  0. 1 13  kg/s. 

Map  flow  ^curacy  is  thought  by  Boeing  engineers  to  be  in  the  order  of  0.05%,  hence  compatible 
with  the  0.1%  goal. 

2.1.2  NEL  facility 

Although  seemingly  less  accurate  than  CEESI  facility,  NEL  gravimetric  system  is  presented  hereafter 
in  order  to  show  some  details  of  the  procedure  which  should  be  followed. 


2. 1.2.1  Test  rig  and  instrumentation 

The  primary  gravimetric  gas  flow  standard  which  is  used  for  the  calibration  tests  is  shown  in  line 
diagram  form  in  fig.8. 

Air  IS  supplied  at  a  pressure  of  210  bar  by  two  reciprocating  compressors  which  have  a  combined  free 
air  deliveiy  of  some  0.06  m^/s.  After  leaving  the  compressors  the  air  passes  through  a  purification 
plant  m  which  ml,  water  vapour  and  solid  particles  are  removed.  The  purification  plant  is  designed  to 
provide  air  with  a  dew-point  not  higher  than  -40°C  and  an  oil  content  not  exceeding  five  parts  per 

million.  A  particle  filter  removes  all  solid  matters  with  size  ratings  greater  than  5  qm. 

This  jprified  air  is  used  to  charge  the  air  storage  vessel  and  the  control  loop  and  in  the  tests  reported 
here  the  air  storage  vessel  is  charged  at  pressures  up  to  76  bar  with  the  control  loop  being  charged  at 
pressures  up  to  65  bar.  During  each  test  run,  air  from  the  control  loop  passes  through  the  pressure 
control  s^ystem  B  via  valve  LI,  which  is  used  to  start  and  stop  the  flow,  to  the  nozzle  under  test 
1  he  air  flowing  from  the  loop  is  replaced  by  air  from  the  air  storage  vessel  which  flows  into  the  loop 
through  pressure  control  system  "A".  Both  of  these  pressure  control  systems  utilise  dome-loaded 
valves  which  are  set  to  give  selected  downstream  pressures. 
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After  passing  through  the  nozzle  under  test  the  air  enters  the  diverter  unit.  In  the  diverter  unit  two 
50  mm  ball  valves  are  located  as  shown  in  fig.8.  When  the  diverter  unit  is  operated  one  valve  closes 
and  the  other  opens  so  that  the  air  flow  can  be  directed  either  through  the  disconnecting  fixture  to  the 
spherical  weighing  vessel  or  through  the  test  and  bypass  lines  to  atmosphere.  The  weighing  vessel 
disconnect  fixture  consists  of  valves  and  a  scaling  system  so  arranged  that  the  charged  weighing 
vessel  can  be  completely  disengaged  froin  the  system  before  being  weighed. 

The  weighing  vessel  is  a  hollow  steel  sphere  1.5  m  in  diameter.  It  is  capable  of  withstanding 
pressures  up  to  300  bar  at  temperatures  in  the  range  -20  to  +50°C.  To  weigh  the  vessel  and  its 
contents  a  sensitive  platform  scale  is  used.  The  scale  is  back-balanced  to  allow  for  the  tare  weight  of 
the  spherical  vessel.  The  weights  are  read  on  a  steelyard  indicator  and  the  sensitivity  of  the  system  is 
such  that  changes  of  mass  of  10  g  can  be  detected  in  diverted  masses  of  typically  40  to  60  kg. 

The  instrumentation  used  during  the  tests  is  as  follows  : 

Nozzle  under  test  ■  ,  • 

The  probe  containing  the  nozzle  upstream  pressure  and  temperature  sensors  is  designed  to  measure 
the  stagnation  values  of  these  parameters.  The  following  accuracies  are  required  to  meet  the  0.1%  goal 
in  mass  flow  measurement  (see  also  chapter  2.4): 

-  upstream  pressure  (PO):  accuracy  better  than  ±0.02  per  cent  of  reading. 

-  upstream  temperature  (TO):  accuracy  better  than  ±0.25°C. 

IJnweighed  volume  ,  r.  f 

Pressure  :  Wallace  and  Tieman  precision  pressure  gauge  of  range  0-35  bar  or  Barnet  test  gauge  ot 

range  0-80  bar.  Accuracy  :  ±0.3  per  cent  of  reading. 

Temperature  :  Rosemount  platinum  resistance  thermometer  with  ASL  type  F25  measurement  system. 
Accuracy  :  ±0.25°C. 

Weighbridge  :  Avery  type  4206  ABA  of  range  0-500  kg. 

Accuracy  :  ±0.02  kg  of  reading. 

Barometric  pressure  :  Model  145,  Texas  Instruments  Ltd,  precision  pressure  gauge  with  type  1 
capsule  of  range  0-1080  mbar  abs.  installed.  Accuracy  :  ±0.05  per  cent  of  reading. 

Ambient  temperature  :  Mercury-in-glass  thermometer.  Accuracy  :  ±1°C. 


2. 1.2. 2  Test  procedure 

The  weighing  vessel  is  lowered  from  the  disconnect  fixture  onto  the  weighbridge  and  its  "empty"  (i.e. 
with  the  air  inside  at  ambient  pressure)  weight  measured.  This  reading  is  recorded  and  the  temperature 
of  the  air  surrounding  the  weighing  vessel  is  also  noted.  The  vessel  is  then  raised  into  its  operating 
position  in  the  disconnect  fixture,  the  securing  lock  engaged,  valves  HI  and  H2  opened  and  vent 
valve  H3  closed. 

Valve  LI  is  opened  very  gradually  to  initiate  flow  through  the  nozzle  under  test,  valves  Dl,  L2,  L3, 
Bl,  B2,  Cl,  C2,  C3  and  C4  already  having  been  opened. 

The  pressure  upstream  of  the  nozzle  under  test  is  then  set  by  adjusting  the  regulating  valves  of 
pressure  control  system  "B".  When  conditions  have  stabilised  the  diverter  mechanism  is  actuated 
switching  the  flow  into  the  weighing  vessel  and  starting  the  timer.  At  the  instant  of  diversion  the 
pressure  and  temperature  of  the  air  in  the  unweighed  volume  are  recorded.  During  the  diversion 
period  the  pressure  and  temperature  upstream  of  the  nozzle  under  test  are  noted  at  regular  intervals. 
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When  the  pressure  in  the  weighing  vessel  reaches  approximately  80  per  cent  of  the  nozzle  upstream 
pressure  (or  40  bar,  whichever  is  the  lower)  the  diverter  is  again  operated  and  at  the  moment  of 
the  pressure  and  temperature  of  the  air  in  the  unweighed  volume  recorded.  Valves  LI  HI 
is  opened,  the  lock  is  disengaged  and  the  vessel  then  lowered  on  to  the 
weighbndge  to  obtain  its  weight  after  the  test.  This  weight,  the  diversion  time  and  the  temperature  of 
t  e  air  surrounding  the  vessel  are  noted.  During  the  test  run  the  barometric  pressure  is  also  recorded. 

The  vessel  is  raised  once  again  into  the  disconnect  fixture,  the  lock  is  engaged,  valves  HI  and  H2  are 
opened  and  then  valve  H3  is  opened  to  exhaust  the  air  from  the  vessel  in  preparation  for  the  next  test 


The  above  procedure  is  repeated  for  each  test  point  carried  out. 
2. 1.2. 3  Test  results 


The  reference  mass  flowrate  through  the  nozzle  under  test  is  obtained  from  the  gravimetric  system  and 
IS  given  by  iWt  i.e.  the  ratio  of  the  mass  (m)  introduced  in  the  vessel  to  the  duration  (t)  of  the 
diversion  penod.  ^  ^ 


Minor  corrections  are  made  (e.g.  atmospheric  buoyancy  etc. . .)  to  get  m  with  the  best 
the  weighing  machine. 


precision  from 


Nozzles  results  are  given  with  references  to  the 
NEL  states  that  the  overall  accuracy  of  the  process  is  0.25%,  and  will 
0.15%  or  better. 


nozzle  Reynolds  number, 
very  soon  be  improved  to  reach 


22  Calibration  of  massflowmeters  by  using  transfer-standard  nozzles 


2.2.1  Types  of  airflow  transfer  standards. 

Most  wind  tunnel  facilities  use  sonic  venturis  to  measure  air  mass  flow.  The  geometry  of  the  venturi 
is  usually  close  to  that  of  the  Smith-Matz  circular  arc  venturi  described  in  ref .4.  This  geometrv  has  the 
advantage  of  a  very  thin  boundary  layer,  so  that  the  transition  process  between  the  laminar  and 
turbulent  flow  regimes  introduces  a  minimum  of  uncertainty  in  the  discharge  coefficient.  However 
based  on  theoretical  predictions  (ref.  5  &  6),  errors  of  0.3  %  or  greater  in  discharge  coefficient  should 
be  expected.  Careful  analysis  of  the  flow  can  reduce  this  uncertainty  (see  chapter  3).  Nevertheless 
calibrations  traceable  to  primary  standards  are  required  if  much  greater  accuracy  is  to  be  achieved. 

Priraiary  standards  are  designed  to  provide  the  high  accuracy  required  for  calibration  of  transfer 
standards.  As  seen  above,  their  mass  flow  range  is  not  large  enough  to  calibrate  directly  most  of  the 
flowmeters  used  in  wind  tunnel  tests  and  the  methods  are  too  time  consuming  and  expensive  for  the 
routine  calibration  of  flowmeters.  Therefore,  any  facility  requiring  such  accuracy  and  traceability  must 
provide  and  have  calibrated  a  transfer  standard  suitable  for  calibration  of  the  flowmeters  that  are 
normally  used  for  testing.  Calibration  of  the  transfer  standard  must  include  the  adequate  accuracy  and 
U-aceabihty  of  pressure  and  temperature  instrumentation,  which  must  also  be  used  when  calibrating 
ow  meters  at  the  user  facility.  Failure  to  do  so  will  result  in  airflow  measurements  which  have  lost 
the  accuracy  and  traceability  to  the  primary  standard. 

The  ttansfer  standard  must  be  designed  to  meet  the  needs  of  the  user  facility,  as  well  as  the  capabilities 
and  liiratations  of  the  pnmary  standards  laboratory.  The  transfer  standard  can  be  a  single  critical 
venturi,  a  multiple  cntical  venturi  (MCV)  with  binary  sized  venturis  such  as  fig.9,  or  a  multiple  critical 
uniform  siz^  ventuns  such  as  fig.  1.  Each  of  these  transfer  standards  must  include  an 
effective  upstrearn  flow  distnbutor  so  that  flow  distortions  cannot  affect  the  airflow  measurements. 
Significant  deviations  from  the  Smith-Matz  venturi  configuration,  such  as  shown  in  fig.2  perform 
equally  as  well  as  evidenced  in  the  data  of  reference  7.  &  ’  F 
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2.2.2  Recommendations  for  calibration  of  working  venturis. 

Calibration  of  a  working  venturi  with  a  transfer  standard  can  be  accomplished  with  a  minimum  loss  of 
accuracy  using  the  following  procedures  : 

2.2.2. 1 

Steady  state  conditions  must  be  achieved  at  each  test  condition  for  both  the  transfer  standard  and  the 
working  venturi  before  data  are  recorded.  Temperatures  will  take  much  longer  time  than  pressures  to 
stabilize  due  to  the  Joule  -  Thomson  temperature  effects  caused  by  the  pressure  drop  of  source 
pressure  control  as  well  as  the  pressure  drop  between  the  two  sonic  devices.  For  instance,  a  50  bar 
pressure  drop  between  the  two  choked  nozzles  can  result  in  a  10  degrees  decrease  of  temperature. 
Therefore  many  minutes  can  be  necessary  to  stabilize  temperatures  at  the  beginning  of  a  test  sequence 
or  run.  Heat  exchangers  upstream  of  each  sonic  device  could  eliminate  this  problem,  and  improve  data 
accuracy. 

2. 2. 2. 2 

Pressure  measurements  for  both  the  transfer  standard  and  the  working  venturi  can  be  measured  with 
an  accuracy  of  ±  0.01  %  to  ±  0.02  %  at  all  conditions.  A  method  which  can  achieve  this  accuracy 
utilizes  dead  weight  testers,  small  range  differential  transducers,  and  solenoid  valves  as  shown  in 
figure  10.  Transducer  zero  can  be  recorded  at  each  pressure  level  by  switching  the  solenoid  valves  so 
that  both  sides  of  the  transducer  are  exposed  to  the  dead  weight  tester,  then  switching  the  measuring 
side  to  read  the  small  pressure  differences  from  the  dead  weight  testers  to  the  transfer  standard  or  the 
working  venturi  for  the  data  record. 

2. 2. 2. 3 

Temperatures  for  both  the  standard  and  the  working  venturi  must  be  verified  to  an  accuracy  of 
+  0.1°C  or  better  at  ambient  temperature,  as  described  under  the  section  on  instrumentation  calibration 
and  checks  (§2.4.1). 

2.2.3  Boeing  transfer  standard  methodology 

Reference  7  describes  the  Boeing  airflow  transfer  standard  with  accuracy  of  0.05  %  or  better  (figure 
1).  This  standard  is  traceable  to  the  Colorado  Engineering  Experiment  Station,  Inc.,  (CEESI) 
Gravimetric  Mass-Time  System  (  ref.8),  which  is  directly  traceable  to  the  U.S. National  Bureau  of 
Standards. 

The  mass-time  system  is  limited  to  a  maximum  flow  rate  of  0.1 13  kg/s  (0.25  Ibm/sec).  Therefore,  the 
transfer  standard  contains  162  identical  venturis  (fig  2),  each  calibrated  over  the  range  from  0.0182  to 
0.113  kg/s  (0.04  to  0.25  Ibm/sec),  permitting  the  transfer  standard  to  provide  calibrations  from 
0.0182  to  18  kg/s  (0.04  to  40  Ibm/sec).  This  transfer  standard  has  the  added  advantage  that  the 
venturis  operate  only  in  the  laminar  flow  regime,  so  that  when  used  to  calibrate  a  working  venturi,  the 
uncertainty  encountered  during  transition  of  the  working  venturi  can  be  accurately  measured.  The 
CEESI  calibration  of  the  Boeing  transfer  standard  is  described  in  ref.9. 

Essential  to  the  development  and  calibration  of  the  Boeing  transfer  standard  was  the  venturi  screening 
facility  (fig.3).  This  facility  consisted  of  two  venturi  stations  in  series,  totally  submerged  in  a 
circulating  water  bath  maintained  at  26.6  °C  (80  °  F).  This  was  to  ensure  that  the  plumbing,  the 
venturis,  and  the  air  passing  through  the  system  were  all  at  the  same  temperature.  The  upstream 
venturi  was  a  common  reference  venturi  with  a  nominal  throat  diameter  of  0,0067  m  (0.263  in.)  and 
the  downstream  venturi  station  permitted  each  of  the  163  venturis  to  be  installed  in  turn.  The  design 
throat  diameter  was  0.00793  m  (0.313  in.)  for  each  of  the  163  venturis. 

Testing  at  Boeing  with  the  venturi  screening  facility  verified  a  high  level  of  repeatability  (  fig.4),  and 
the  ability  to  accurately  measure  differences  between  the  163  venturis.  Because  the  venturis  were 
manufactured  to  be  as  identical  as  possible,  it  was  not  surprising  that  they  had  virtually  identical 
Reynolds  Number  characteristics,  and  repeatable  but  different  levels  of  discharge  coefficient.  These 
differences  were  found  to  be  caused  by  the  inability  to  accurately  measure  the  venturi  throat  diameter 
(area).  Therefore,  the  throat  areas  were  adjusted  to  reflect  effective  area  rather  than  measured  area. 


which  resulted  in  a  common  Cd  vs  Rd  curve  for  all  163  venturis  (fig  5).  This  process  is  explained  in 
detail  in  ref.7.  This  ability  to  accurately  measure  the  venturi  differences,  eliminated  the  need  to 
calibrate  all  163  venturis  by  gravimetry.  Such  a  calibration  would  have  taken  an  estimated  4  years. 

The  calibrations  at  CEESI  were  done  with  the  Boeing  venturi  screening  facility  installed  directly 
downstream  of  the  CEESI  mass-time  system  (fig.6).  Fourteen  venturis  were  selected  for  calibration 
starting  with  S/N  1  and  ending  with  S/N  160,  with  the  selections  somewhat  uniformly  spaced  to 
cover  the  entire  manufacturing  sequence.  S/N  1  venturi  was  calibrated  at  the  beginning  and  end  of  the 
calibrations,  and  twice  more  between  the  other  venturi  calibrations.  Of  course,  the  common  reference 
venturi  was  calibrated  17  times,  since  it  was  in  series  with  each  of  the  other  venturis.  Using  the 
effective  throat  areas  for  each  venturi  which  were  evaluated  earlier  at  Boeing,  the  CEESI  gravimetric 
calibrations  verified  the  common  Cd  vs  Rd  characteristic  as  well  as  the  absolute  levels  of  Cd  (fig.7) 
for  all  163  venturis. 

An  uncertainty  analysis  considering  all  elements  of  the  calibration  setup  is  presented  in  ref.7.  This 
analysis  estimates  the  uncertainty  of  discharge  coefficient  to  be  ±  0.07  %  for  an  individual  venturi, 
and  ±  0.05  %  for  2  or  more  venturis  in  parallel. 

S/N  1  venturi  has  been  preserved  as  a  permanent  control  venturi  for  future  calibration  work.  It  is  not 
used  as  part  of  the  transfer  standard.  Should  any  of  the  other  162  venturis  become  damaged  or 
deteriorated,  the  venturi  screening  facility  could  be  reactivated  using  S/N  1  venturi  to  verify  the 
integrity  of  the  common  reference  venturi,  which  in  turn  could  evaluate  the  Cd  vs  Rd  characteristic 
and  effective  throat  area  for  replacement  venturis. 


23  Nozzle  thrust  measuring  benches. 


Thmst  measuring  benches  are  required  to  enable  the  thrust  and  discharge  characteristics  of  nozzle 
systems  to  be  determined  in  isolation.  A  number  of  facilities  exist  world- wide  to  either  support  nozzle 
system  design  and/or  for  calibration  of  nacelle  exhaust  nozzle  configurations  for  use  in  wind  tunnel 
testing. 

The  basic  requirements  for  thrust  benches  are  common  to  most  applications.  A  means  is  required 
whereby  the  thrust  of  the  individual  nozzles  (primary  and  fan)  can  be  derived  from  the  measured 
overall  thmst  of  the  system  and  related  to  known  upstream  conditions.  The  nozzle  system  must  be 
placed  in  a  quiescent  exhaust  environment  but  under  conditions  that  match  the  required  model  external 
pressures.  The  test  technique  necessitates  careful  bookkeeping  of  all  the  mass  flows  and  forces 
involved  in  the  system  if  the  required  accuracies  are  to  be  achieved. 

In  general  the  thmst  bench  will  comprise  a  live  (metric)  model  support  frame,  capable  of  measuring  as 
a  minimum  the  force  along  the  nozzle  thrust  axis,  surrounded  by  plenum  shells  and  a  means  to 
provide  the  required  inlet  and/or  exhaust  pressure  environments.  Careful  attention  has  to  be  given 
during  facility  design  to  minimise  the  tare  loads  generated  on  the  live  frame  by  the  transfer  of  air 
across  the  force  balance  system  as  well  as  to  minimise  the  pressure  area  terms  acting  on  it. 

Sonic  venturis  are  typically  used  to  measure  the  air  mass  flows  entering  into  and  exiting  from  the 
system:  these  require  high  precision  and  ideally  should  be  traceable  to  recognised  standards. 

Procedures  are  required  to  routinely  check  correct  operation  of  all  the  facility  instmmentation  (force, 
pressure  and  temperature)  and  to  ensure  that  consistent  calibration  practices  are  followed. 

Nozzles  with  'known'  characteristics  (thrust  and  discharge)  should  be  used  on  a  regular  basis  to 
provide  overall  facility  calibrations  and  to  ensure  continued  satisfactory  operation. 

For  Turbo-Powered-Simulator  tests  dryness  for  primary  and  secondary  flows  must  be  strictly 
controlled  in  order  to  avoid  any  ice  formation.  For  all  tests  air  should  be  dry  enough  to  make  sure  that 
no  risk  of  local  condensation  exists  at  slightly  supersonic  speeds  in  the  nozzles. 
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Different  types  of  thrust  benches  routinely  used  at  ONERA,  DLR,  ARA  ,  NLR  and  Boeing  are 
described  in  appendix  A  to  E. 


2.4  General  precautions  for  all  tests 


2.4.1  Instrumentation  calibration  and  recommended  checks  for  all  tests 

Nothing  can  be  left  to  chance  if  the  goal  of  0.1%  accuracy  of  thrust  measurement  is  to  be  achieved  or 
nearly  achieved.  The  only  item  that  is  not  calibrated  or  checked  in  the  test  set-up  will  probably  be  the 
one  that  does  not  work  as  advertised.  There  is  nothing  so  wasteful  as  running  a  test  with 
undiscovered  and  uncorrected  errors  or  malfunctions.  Calibrations  and  checks  should  be  done  prior  to 
every  test,  and  repeated  at  intervals  during  a  test  to  develop  confidence  in  the  reliability  and  stability 
of  the  data.  The  intervals  can  be  :  every  shift ;  every  day  ;  or  anytime  a  malfunction  is  suspected. 

Pressure  transducers : 

Must  be  isolated  from  acoustic  noise,  vibration,  thermal  input. 

Must  be  calibrated  with  a  traceable  standard-dead  weight  tester. 

Must  exhibit  small  hysteresis  characteristics. 

Curve  fit  of  data  must  be  within  ±  0.02%.  (1) 

Check  calib's  each  shift  must  repeat  within  ±  0.02%.  (1) 

Force  Measurement  Systems : 

Precision,  traceable  weights  used  for  all  loads. 

All  components  of  the  balance  must  be  calibrated. 

Horizontal  forces  use  a  cable  and  pulley  or  knives  for  loading. 

Axial  force  pulley  must  be  at  least  12"  dia.  low  friction. 

Curve  fit  of  axial  force  data  must  be  within  ±0.02%.  (1) 

Check  calib's  each  shift,  axial  force  must  repeat  within  ±  0.02%  .  (1) 

Other  5  components  can  be  less  accurate. 

Temperature  Sensors : 

When  placed  in  a  suitable  heat  sink  with  a  traceable  standard  thermometer,  should  agree  with  the 
standard  within  ±  0.1°C  or  better  after  stabilization.  (2)  or  (3) 

Traceable  Pressure  Standards  ;  Dead  weight  testers. 

Make  sure  Metrology  Lab.  certifications  meet  your  needs. 

Do  not  accept  certifications  at  face  value. 

Understand  Metrology  Lab.  methods  of  data  analysis. 

(1)  ±  0.02%  of  data  values  in  the  upper  80%  of  range. 

(2)  For  removable  sensors,  an  aluminium  block  makes  an  excellent  heat  sink,  with  holes  to  insert 
sensors  and  a  standard  thermometer.  When  placed  on  a  styrofoam  block  away  from  draughts  and  heat 
sources,  it  will  stabilize  to  room  temperature  in  15  to  20  minutes. 

(3)  For  sensors  permanently  mounted  in  a  blowing  nacelle  or  Turbo-Powered-Simulator,  plug  the 
inlet/exhaust  with  soft  cloths,  carefully  insert  two  or  more  standard  thermometers  past/through  the 
cloths  so  the  sensing  bulbs  are  in  contact  with  structure  near  the  sensors.  Wrap  the  entire  nacelle  with 
a  blanket  so  the  only  source  of  thermal  input  is  through  the  strut.  Protect  from  draughts  and  heat 
sources.  Allow  to  stabilize. 
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2.4.2  Determination  of  Tare  Loads 

2.4.2. 1  General 

For  tests  on  powered  models  (isolated  or  installed)  compressed  air  is  often  used  as  power  source.  The 
compressed  air  has  to  be  supplied  to  the  powered  model  through  an  air  duct  which  has  to  pass  the 
balance  system  with  only  minor  interaction.  In  the  calibration  of  the  (internal  or  external)  balance 
system,  only  dead  weight  effects  and  elasticity  effects  of  the  inoperative  ( decompressed)  air  duct  have 
been  taken  into  account  implicitly.  When  compressed  air  and  possibly  heated  air  flows  through  the  air 
duct  towards  the  model  in  the  test  section,  the  duct  may  exert  residual  forces  and  moments  on  the 
balance  structure,  that  will  be  interpreted  by  the  balance  system  as  forces  and  moments  in  the  Balance 
Centre  BC.  Three  separate  effects  in  the  air  duct  can  be  distinguished,  although  they  occur 
simultaneously  when  compressed  air  is  supplied  to  the  powered  model : 

2.4.2. 1.1  Pressure 

The  pressure  inside  the  air  duct  deforms  the  structure  in  such  a  way,  that  an  extra  force  may  act  on  the 
balance  system.  This  effect  depends  on  the  magnitude  of  the  pressure  inside  the  duct.  A  part  of  this 
pressure  effect  is  also  due  to  the  weight  of  the  compressed  air.  The  type  of  uncoupling  bellows  and 
the  overall  layout  of  the  rig  can  considerably  reduce  any  of  these  effects. 

2.4.2. 1.2  Temperature 

When  the  air  duct  is  heated,  the  duct  expands  and  consequently  an  extra  force  may  act  on  the  balance 
system.  This  effect  depends  mainly  on  the  temperature  diference  between  air  duct  and  balance 
structure. 

2.4.2. 1.3  Momentum  tare 

In  order  to  avoid  any  momentum  effect  of  the  air  at  the  entrance  of  the  metric  part  of  the  balance,  the 
momentum  of  the  air  coming  from  the  earth-fixed  ducts  and  entering  the  ducts  of  the  metric  part  of  the 
balance  must  be  strictly  perpendicular  to  the  axis  of  the  nozzle. 

The  non  fulfilment  of  this  conditions,  e.g.  if  the  velocity  disk  in  the  air  duct  is  possibly  skew,  gives 
rise  to  a  momentum  effect.  This  effect  depending  on  both  pressure  and  velocity  levels  at  the  entrance 
of  the  metric  part  of  the  balance  must  be  carefully  evaluated. 

The  air  duct  pressure  and  temperature  effects  on  the  balance  system  can  easily  be  determined 
statically.  Proper  evaluation  of  the  momentum  tare  loads  is  on  the  other  hand  more  cumbersome  but 
can  be  accomplished  with  reference  nozzles.  Two  different  methods  (nozzles)  will  be  discussed 
beneath.  An  alternative,  using  standard  nozzles,  is  presented  in  ref  [28] . 


2. 4. 2. 2  Zero-thrust  body 

Evduation  of  momentum  tare  can  be  accomplished  with  a  zero  -  thrust  -  body,  as  shown  in  fig.  1 1. 
This  is  a  simple  device,  but  it  must  be  accurately  constmcted  so  that  the  nozzle  centrelines  are 
coincident  with  each  other,  and  perpendicular  to  the  main  body  centreline.  The  bolt  pattern  at  each  end 
of  the  main  body  should  be  identical,  so  that  the  body  can  be  reversed  end  -  to  -  end,  or  rotated  about 
its  centreline  in  90  °  increments.  A  thick,  subsonic,  perforated  plate  is  needed  at  the  entrance  to  the 
tube  supporting  each  nozzle  to  align  and  distribute  the  flow.  A  series  of  pairs  of  nozzles  is  required  so 
that  4  or  5  sizes  can  be  used  to  evaluate  each  size  airbridge.  The  nozzles  should  be  sonic  at  all 
conditions,  which  keeps  the  air  duct  at  a  constant  velocity  whatever  be  the  pressure  level,  for  each 
nozzle  size.  A  pressure  and  temperature  measurement  is  required  at  the  end  of  the  body  as  shown. 

Once  the  nozzle  pairs  have  been  calibrated  by  a  sonic  venturi,  the  nozzles  can  be  used  to  measure  the 
mass-flow  rate  for  momentum  tare  evaluation,  using  the  pressure  and  temperature  in  the  zero-thrust 
body. 
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Each  of  the  nozzle  pairs  may  not  be  identical  in  size,  and  as  a  result,  may  not  cancel  or  produce  zero 
thrust.  Therefore,  it  is  recommended  that  when  the  nozzles  are  oriented  in  the  plane  of  balance  side 
force,  that  both  side  force  and  yaw  moment  loads  are  ignored.  Similarly,  when  the  nozzles  are 
oriented  in  the  plane  of  normal  force  (vertical  plane),  that  both  normal  force  and  pitching  moment 
loads  are  ignored. 


For  each  powered  test  installation,  a  pressure  and  temperature  measurement  must  be  installed  at  the 
discharge  end  of  the  air  duct  system  to  relate  to  the  pressure  and  momentum  tare  equations.  These 
equations  are  shown  in  figure  1 1 .  Total  mass  flow  for  the  air  duct  can  be  obtained  from  the  sonic 
venturi/venturis  used  for  the  test. 


2. 4. 2. 3  Blown  nozzle  in  the  wind  tunnel 

An  evaluation  of  momentum  tare  in  the  wind  tunnel  can  also  be  accomplished  with  a  blown  nozzle. 
The  blown  nozzle  should  deliver  thrust  along  the  nozzle  centreline.  The  blown  nozzle  can  be  used 
either  in  an  isolated  or  installed  test  set-up.  The  principle  is  explained  below  for  an  isolated  blown 
nozzle  which  is  mounted  to  a  strut  (air  duct)  for  air  supply  and  connection  to  the  external  balance 
system. 

The  strut-nozzle  interface  should  be  designed  in  such  a  way  that  the  nozzle  can  be  mounted  to  the  strut 
in  "up-wind"  and  "down-wind"  blowing  direction.  In  this  way  the  thrust  exerts  the  external  balance 
during  one  test  in  positive  and  during  the  other  test  in  negative  direction. 

The  isolated  blown  nozzle  should  be  tested  at  wind-off  conditions  to  determine  the  internal  thrust  and 
mass  flows.  The  exhaust  flow  however  can  entrain  flow  completely  around  the  tunnel  circuit  leading 
to  external  aerodynamic  forces  on  the  metric  parts  (  strut,  nozzle,  etc.).  This  can  be  suppressed  as 
follows : 

a)  Main  jet  flow  in  tunnel  direction  (down  -  wind  blowing,  see  fig.  12  A). 

The  tunnel  test  section  is  closed  by  a  barrier  ahead  of  the  test  section,  surplus  air  should  be  vented  to 
the  environment  at  the  settling  chamber. 

b)  Main  jet  flow  in  opposite  tunnel  direction  (up  -  wind  blowing,  see  fig.  12  B). 

The  tunnel  test  section  is  closed  by  a  barrier  behind  the  test  section,  surplus  air  should  vented  to  the 
environment  at  the  diffusor  section  (not  shown  in  fig.  12  B). 

A  second  non  -  metric  barrier  might  be  required  just  upstream  of  the  nozzle  exit,  to  confine  secondary 
recirculation  flow  downstream  of  this  barrier  and  separate  from  the  metric  parts. 

The  blown  nozzle  should  be  supplied  with  pressure  and  temperature  rakes.  During  the  test 
compressed  air  is  supplied  to  the  nozzle  up  to  the  operational  levels  of  pressure  and  flow.  The  mass 
flow  can  be  obtained  from  the  sonic  venturi  mounted  in  the  drive  air  supply  line.  The  balance  readings 
should  have  been  properly  corrected  for  the  static  pressure  and  temperature  effects  (as  described 
above)  and  for  the  measured  strut  deflections  (  change  of  engine  centre  N  with  respect  to  balance 
centre  BC  and  rotation  of  nozzle  axis  respectively),  which  are  introduced  by  forces,  moments, 
temperature  effects  etc. 

The  blown  nozzle  should  be  calibrated  in  down-wind  and  up-wind  direction  and  the  derived 
components  of  the  engine  velocity  coefficient  vector  are  determined  as  a  function  of  nozzle  pressure 
ratio  NPR. 

The  momentum  tare  loads  are  derived  from: 

1)  The  differences  in  the  magnitude  of  the  velocity  vector  components  between  both  test  set  -  ups  (in 
this  way  the  force  and  moment  in  Z  -  direction  is  not  derived). 
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2)  The  magnitudes  of  the  velocity  vector  components,  with  exception  of  X-  direction.  Theoretically  all 
these  thrust  components  should  be  zero. 

One  should  be  aware  that  only  sensible  momentum  tare  loads  are  derived  when  they  are  larger  than 
might  be  expected  from  the  balance  accuracy  and  nozzle  misalignment  errors. 

Some  explanations  concerning  this  technique  of  momentum  tare  evaluation  with  down-wind  and  up¬ 
wind  blowing  nozzles  are  given  in  figure  13. 
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3.  FLOW  ANALYSIS 


There  are  several  good  reasons  for  paying  attention  to  the  detailed  behaviour  of  the  flow  in  the  two 
types  of  nozzles  which  are  currently  in  use  for  mass  flow  and  thrust  measurements  in  aerodynamic 
facilities 

a-  it  may  be  difficult  to  justify  the  time  and  expense  that  Boeing  expended  in  developing  its 
Transfer  Standard  tools  and  methodology.  However,  the  sanae  repeatability  and  consistency  of  mass 
flow  measurement  can  be  achieved  by  calibrating  each  working  venturi  in  series  with  a  theoretically 
defined  Standard  venturi  even  though  absolute  accuracy  may  be  slightly  in  error.  Cross  checking  data 
with  other  facilities  can  determine  if  a  systematic  error  remains  in  the  absolute  level  of  mass  flow. 
This  implies  that  accuracy  of  theoretical  predictions  is  much  better  than  the  above  quoted  0.3  %. 

b-  assumptions  of  how  toroidal  throat  nozzles  do  work,  i.e.  which  boundary  layer  thickness 
should  be  accounted  for  at  the  throat,  were  slightly  different  in  the  different  establishments,  generating 
differences  in  the  order  of  0.2%  on  calculated  discharge  coefficients. 

c-  most  members  of  the  group  had  had  bad  experiences  in  using  cylindrical  throat  nozzles 
(ASME  LRN)  as  references  and  they  could  not  state  precisely  why;  see  for  example  fig.  32, 33,  34. 
Hence  the  Working  Group  decided  that  it  was  worth  asking  specialists  to  proceed  with  a  few 
calculations  and  physical  discussions  in  order  to  clarify  above  items  (b)  and  (c)  and  to  make  a  clear 
statement  on  what  can  be  achieved  with  (a). 


3.1  Cylindrical  throat  nozzles:  the  ASME  Long  Radius  Nozzle 


3.1.1  General  considerations  on  ASME  Long  Radius  Nozzle 

The  ASME  Long  Radius  Nozzle  (ASME  LRN)  is  composed  of  a  convergent  inlet  section  described 
by  a  quarter  ellipse,  and  a  circular  cylinder  throat  section.  The  connecting  plane  between  the  two 
sections  is  at  their  point  of  tangency.  There  are  two  designations  for  the  ASME  LRN  based  on  internal 
geometric  shape.  The  geometric  shape  is  defined  by  the  ratio  of  the  nozzle  throat  diameter  (d)  to  the 
upstream  internal  pipe  diameter  (D).  This  ratio  is  denoted  as  "p",  where  P  =  d/D.  Typically  p  falls 
within  the  range  of  0.2  <  P  <  0.8.  When  p  >  0.45  the  nozzle  is  referred  to  as  a  "high  p"  nozzle  and 
when  p  <  0,5  it  is  referred  to  as  a  "low  p"  nozzle.  The  range  between  0.45  and  0.50  is  common  to 
both  types  of  nozzles.  The  relationships  between  p  and  the  various  geometric  parameters  describing 
the  convergent  inlet  section,  the  cylindrical  throat  section  and  the  internal  pipe  and  throat  diameter  are 
shown  in  figure  14.  From  a  practitioners  perspective  the  "high  p"  nozzle  minimizes  flow  restriction 
within  the  piping  system  while  the  use  of  a  "low  p"  nozzle  causes  a  greater  flow  restriction  for  a  given 
mass  flow  and  thus  a  correspondingly  higher  pressure  differential  which  yields  greater  accuracy. 

The  ASME  Long  Radius  Nozzle  was  designed  for  and  is  traditionally  operated  with  subsonic  flow  in 
the  nozzle  throat.  However,  the  nozzle  can  also  be  operated  in  a  critical  or  choked  flow  condition,  in 
which  case  the  massflow  rate  of  a  given  fluid  through  a  given  nozzle  can  in  principle  be  easily 
calculated  from  simple  one-dimensional  theory  as  it  is  only  a  function  of  upstream  total  temperature 
and  pressure.  This  theoretical  massflow  rate  however  is  always  different  from  the  actual  one  due  to  ; 

1-  real  gas  effects; 

2-  growth  and  development  of  the  boundary  layer  within  the  nozzle; 

3-  sonic  line  curvature. 

The  coefficient  of  discharge  (Cd)  is  defined  as  the  ratio  of  the  actual  to  the  theoretical  massflow  rate 
and  is  expected  to  account  for  all  discrepancies  between  theory  and  the  actual  flow  conditions. 
Examples  of  discrepancies  are  well  documented  by  B.T.  Amberg  [14]. 

Very  few  references  pertaining  to  the  coefficients  of  discharge  and  the  operation  of  ASME  LRN  at 
high  subsonic  and  critical  or  choked  flow  conditions  are  available  in  the  open  literature.  Those  that  are 
available  indicate  that  the  uncertainty  in  the  values  of  Cd  was  significant  when  the  nozzle  operates 
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With  high  subsonic  or  critical  throat  flow.  This  uncertainty  was  primarily  attributed  [17, 18, 19,  20]  to 
boundary  layer  transition  which  occurred  at  throat  Reynolds  numbers  between  0.6E6  to  2E6.  The 
transition  was  promoted  by  the  local  adverse  pressure  gradient  generated  by  the  discontinuity  in  the 
wall-radius  of  curvamre  at  the  juncture  of  the  elliptical  inlet  section  and  the  cylindrical  throat  section. 
The  sensitivity  of  this  boundary  layer  transition  to  boundary  layer  development  in  the  converging  inlet 
section  was  discussed  by  Reimer  [15]  who  found  that  the  effects  of  machining  and  polishing  could 
increase  the  Cd  by  0.0025  for  a  given  nozzle.  Smith  and  Matz  [17]  found  for  critical  flow  the 
uncertainty  in  Cd  resulting  from  transition  was  further  exacerbated  by  the  presence  of  a  local 
supersonic  flow  zone  at  the  juncture  of  the  elliptical  inlet  section  and  the  cylindrical  throat  section. 
This  supersonic  flow  region  was  likely  accompanied  by  downstream  lambda  shocks.  All  the  above 
flow  phenomena  affect  the  level  of  uncertainty  in  Cd  due  to  sensitiveness  of  boundary  layer  and  flow 
development  within  the  nozzle.  The  magnitude  of  this  uncertainty  was  of  the  order  of  ±0.3%  (95% 
confidence  band)  over  the  transition  range.  Uncertainty  was  reduced  to  +  0.15%  downstream  of  the 
transition  zone. 

In  order  to  further  understand  the  physics  of  compressible  air  flow  through  a  low  6  value  ASME 
LRN,  a  2-D  axisymetric  CFD  study  was  undenaken  by  Kubberud  of  CFD  norway.  Both  Euler  and 
viscous  Navier-Stokes  codes  were  employed.  The  Euler  code  (1515  grid  points)  was  run  for  two 
nozzle  geometries,  first  an  ASME  LRN  and  second,  a  modified  ASME  LRN  that  is  with  a  divergent 
exit  section  mated  to  the  downstream  end  of  the  cylindrical  portion  of  the  nozzle.  The  grids  for  the 
two  geometries  are  shown  in  fig.  15. 

The  Navier-Stokes  code  (6161  grid  points,  Baldwin  Lomax  turbulence  model,  Sutherland's  viscosity 
formula,  Red=8  10^)  was  run  using  the  second  geometry. 

Upstream  boundary  conditions  were  as  follows;  total  pressure  :  45  bar;  total  temperature:  293K. 
Figures  16a  and  16b  (Euler  code)  and  16c  (Navier-Stokes  code)  show  computed  iso-Mach  number 
lines  superimposed  on  the  nozzle  outlines.  For  the  two  Euler  calculations,  the  location  of  the  sonic 
line  was  found  to  be  independent  of  nozzle  exit  geometry  as  is  shown  fig.  16a  and  16b.Figure  16c 
shows  the  effect  of  viscosity  (growth  of  the  wall  boundary  layer)  on  the  location  of  the  sonic  line 
which  has  moved  approximatively  one-half  a  nozzle  diameter  downstream  from  the  location  shown  in 
fig.  16a  and  16b.  A  reduction  of  0.8%  in  the  discharge  coefficient  Cd  results  from  the  growth  of  the 
boundary  layer.  The  table  hereafter  shows  the  results. 


)  sonic  line  curvature  effect  only 

) 

)  boundary  layer  and  sonic  line 
)  curvature  effects 

These  effects  also  exist  in  other  nozzles  such  as  short  radius  nozzles  discussed  below.  The  point  is 
that  they  are  about  twice  stronger  with  the  ASME  LRN. 

Accurate  determination  of  Cd  for  a  critical  flow  nozzle  is  only  one  of  the  many  parameters  required  to 
calculate  the  mass  flow,  see  Amberg  [14].  The  uncertainty  in  Cd  of  the  ASME  LRN  exceeds  the 
desired  uncertainty  of  0.1%  in  the  measurement  of  the  mass  flow  rate  required  by  AG ARD  Working 
Group  19.  The  experimental  data  cited  in  the  literature  [15  to  21]  typically  He  0.7%  below  the  ASME 
curve  and  the  strong  experimental  evidence  of  boundary  layer  transition  is  neither  shown  nor 
replicated  by  the  ASME  curve. 

This  uncertainty  makes  the  ASME  Long  Radius  Nozzle  unsuitable  for  use  as  a  sonic  nozzle  for  mass 
flow  measurement  and  thrust  calibration.  The  ASME  Long  Radius  Nozzle  was  not  designed  for 
operation  with  high  subsonic  or  critical  flow  conditions  within  the  nozzle  throat.  The  fact  that  its 
ability  to  measure  mass  flow  rate  with  minimal  uncertainty  is  unsatisfactoiy  reflects  not  on  the  nozzle 
design  but  on  the  inappropriate  operating  condition  of  the  nozzle. 


Flow  condition 

Cn 

Isentropic 
(one  dimensional) 

1.0 

Non  viscous,  Euler 
(Axisymetric) 

0.99916 

Viscous,  Navier-Stokes, 
Red=  8 106 
(Axisymetric) 

0.99088 
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32  Toroidal  throat  nozzle 

3.2.1  Short  radius  nozzles:  descriptions  and  references 

In  wind  tunnels  or  other  facilities  where  accurate  mass  flow  measurements  are  required,  mass  flows 
are  usually  measured  through  sonic  venturis  with  short  radius  nozzle.  This  geometry  has  the 
advantage  of  a  well  known  curvature  effect  and  of  a  very  thin  boundary  layer. 

Smith-Matz  recommends  a  circular  arc  nozzle  (fig.  17),  (ref.4);  the  ratio  of  the  arc  radius  to  the  throat 
diameter  is  1.8175.  ONERA  uses  similar  circular  arc  nozzle  with  a  ratio  of  2  (Appendix  F).  For  the 
Boeing  transfer  standard  nozzle,  Stevens  (ref.7),  uses  a  continuous  curvature  entrance  shape  with 
such  parameters  that  in  the  throat  region,  the  curve  shape  approximates  the  Smith-Matz  circular  arc 
contour  (fig.  17). 

The  centrifugal  forces  created  by  the  turning  of  the  flow  in  the  contraction  section  produces  a  non  one 
dimensional  flow  and  then,  a  non  uniform  pressure  and  velocity  distribution  at  the  throat.  This  effect 
can  be  accurately  predicted  for  circular  arc  nozzles  (ref.  23  and  25). 

The  area  reduction  due  to  the  boundary  layer  can  be  also  accurately  calculated,  but  the  nature  of  the 
boundary  layer  (laminar  or  turbulent)  can  produce  differences  of  0.2  to  0.3%.  The  nature  and 
thickness  of  the  boundary  layer  at  the  throat  are  subject  to  a  discussion  beneath. 

3.2.2  Examples  of  discharge  coefficient  calculations  (calorically  perfect  gaz,  y=  1.4  ) 

Various  numerical  flow  analysis  have  been  undertaken  on  the  ONERA  short  radius  nozzle  described 
in  appendix  F,  by  CFD  norway  at  Trondheim  and  by  the  Aerodynamic  Division  at  ONERA. 

Four  cases  of  calculations  can  be  reported  and  summarized  as  follows: 


CASE 

Throat  diameter  (mm) 

Stagnation  pressure 

(bar) 

Reynolds  number 

Red 

1 

10 

25 

3.8  106 

2 

10 

45 

6.8  106 

3 

20 

25 

7.6  106 

4 

20 

45 

14  106 

The  Reynolds  number  Red  is  calculated  with  the  throat  diameter  and  the  throat  flow  conditions 
(critical  flow): 

ac  d  pc  ^  ^  stagnation  temperature  of  293  K  was  held  constant  in  all  cases. 


Rtd  = 


Pc 


3.2.2. 1  CFD  norway  results 

The  associated  numerical  grid  for  non-viscous  and  viscous  calculations,  i.e.  the  grid  used  for  solving 
the  axisymetric  Euler  and  Navier-Stokes  equations  respectively,  are  shown  in  figure  18.  The 
discharge  coefficient  Cd  has  been  evaluated  with  the  results  tabulated  in  the  following: 


Flow  conditions 

Discharge 
coefficient  CD 

Case  1 

Discharge 
coefficient  Cd 

Case  2 

Discharge 
coefficient  Cd 

Case  3 

Discharge 
coefficient  Cd 

Case  4 

Isentropic  1-dim 

1.0 

1.0 

1.0 

1.0 

Non  viscous 

0.99889 

0.99888 

0.9983 

0.9984 

Viscous,  laminar 

0.99685 

0.99733 

0.9973 

- 

Viscous  tuibiilent 

0  99330 

iill 

0.99472 

Ill* 

0.99500 

J 

Viscous 

relaminarized 

0.99330 

0.99359 

- 

- 

Case  3  was  calculated  for  both  laminar  and  turbulent  flow  (Stock-Haase  turbulence  model),  and  the 
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results  are  presented  figure  19.  As  it  can  be  seen,  a  laminar  separation  bubble  is  visible  in  the  first  part 
of  the  nozzle.  However,  the  influence  on  the  sonic  line  at  the  throat  is  negligible  and  the  turbulent 
flow  case  is  regarded  as  the  most  realistic  one  (Red  =  7.6  10^). 


3. 2. 2. 2  ONERA  results 

The  same  calculations  have  been  performed  independently  by  ONERA.  The  results  are  the  following: 


Flow  conditions 

Discharge 
coefficient  Cd 

Discharge 
coefficient  Cd 

Discharge 
coefficient  Cd 

Discharge 
coefficient  Cd 

Case  1 

Case  2 

Case  3 

Case  4 

Isentropic  1-dim 

1.0 

1.0 

1.0 

1.0 

Non  viscous 

0.9991 

0.9991 

0.99902 

0.99902 

0.99428 

0.99477 

0.99515 

When  making  the  same  physical  assumptions,  calculations  should  give  very  close  results.  In  the 
"simple"  non-viscous  case,  differences  can  however  attain  0.06%  and  in  the  viscous  case  they  can 
slightly  exceed  0.1%  which  is  already  too  much! 


3. 2. 2. 3  Comparison  with  previous  results 

It  is  worth  asking  whether  these  sophisticated  codes  bring  much  improvement  comparing  with 
analytical  calculations  made  long  ago  by  Masure  in  1968  (ref  23)  and  Green  in  1971  (ref  25).  Results 
according  to  Masure  (turbulent  case)  and  Green  ((laminar  case)  are  shown  hereafter,  compared  with 
the  above  results.  They  have  been  derived  respectively  from  ref  23  and  25. 

Nota:  Masure's  charts  extracted  from  ref  30  and  reproduced  figure  17bis,  allow  to  predict  quickly 
short  nozzle  flow  and  thrust  coefficient.  These  charts  will  be  used  for  the  so  called  «  practical 
formulae  »  of  chapter  3.5.  Warning!  Reynolds  number  mentioned  on  these  charts  is  calculated  with 
the  half  throat  diameter  h=d/2,  and  upstream  stagnation  conditions: 


flo 

One  finds:  Red  =  1,34  Reo,h 

Non  viscous 


Sonic  line 
curvature  effect 
only 


Case 

1 

2 

3 

4 

CFD  norway 

0.9989 

0.9989 

0.99833 

0.99833 

ONERA 

0.9991 

0.9991 

0.99902 

0.99902 

Masure 

0.99857 

0.99857 

0.99857 

0.99857 

Green 

0.99844 

0.99844 

0.99844 

0.99844 

Viscous  (turbulent) 


Boundary  layer 
and  sonic  line 
curvature 
effects 


Green's  results  are  not  given  because  they  are  calculated  using  laminar  boundary  layer  at  the  throat. 
One  can  see  that  all  these  results  can  deviate  by  up  to  0.0007  for  the  sonic  line  curvature  effect  and  by 
more  than  0.001  for  the  viscous  (turbulent)  case. 


Red 

3.8  106 

6.8  106 

7.6  106 

14  106 

Isentropic  1-dim 

1 

1 

1 

1 

CFD  norway 

0.9933 

0.99359 

0.99472 

0.99509 

ONERA 

0.99428 

0.99477 

0.99471 

0.99515 

Masure 

0.99475 

0.99510 

0.99516 

0.99548 
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Nevertheless,  assuming  that  this  0.1%  difference  between  calculations  will  be  eliminated  some  day,  it 
was  decided  to  discuss  the  more  important  topic  of  the  boundary  layer  behaviour  in  order  to 
understand  why  it  could  have  been  considered  laminar  by  one  author  and  turbulent  by  another,  the 
result  being  an  uncertainty  of  0.2%  to  0.3%  on  the  predicted  discharge  coefficient. 


3. 2. 2. 4  Boundary  layer  discussion 


It  is  important  for  both  theoretical  predictions  and  experimental  work  to  know  the  nature  of  the 
boundary  layer  i.e.  laminar  or  turbulent  in  the  nozzles: 

-  mass-flow  computations  may  change  by  about  0.2%  if  a  wrong  assumption  is  made  in  theoretical 
models, 

-  selection  of  types  of  nozzles,  machining  and  handling  precautions,  avoidance  of  transitional 
situations  are  directly  affected  by  actual  or  assumed  boundary  layer  conditions. 

Such  a  discussion  is  not  new.  TTie  topic  has  been  adressed  in  the  30's  and  more  in  depth  by  Hall  in 
1959  (ref.27).  However,  since  reputable  scientists  like  Masure  and  Green  have  been  led  on  this  matter 
to  what  at  first  glance  was  appearing  as  contradicting  assumptions,  the  working  group  had  initially 
some  difficulties  to  understand  how  different  approaches  could  result  in  what  seemed  to  be  a 
reasonably  good  and  consistent  set  of  data.  It  was  unfortunately  impossible  to  find  an  unquestionable 
experimental  description  of  the  boundary  layer  behaviour  for  the  considered  family  of  sonic  nozzles 
through  a  large  enough  Reynolds  number  range,  however : 

-  well  below  10^  R.L.Steven's  direct  calibration  with  mass-time  standards  shows  trends  versus 
Reynolds  number  inferring  that  throat  boundary  layers  are  laminar; 

-  well  above  10^  numerous  specialists  following  Masure  have  made  with  some  good  reasons  the 
assumption  that  the  throat  boundary  layer  is  turbulent. 

Suspecting  that  both  might  be  right  at  least  at  their  end  of  the  Reynolds  number  range  and  being  aware 
of  recent  progress  in  the  field  of  laminar  flow  control,  a  few  calculations  of  typical  nozzle  flows  using 
up  to  date  models  have  been  performed  by  boundary  layer  specialists  at  CFD  norway  and  at  ONERA- 
CERT. 

Difficulties  start  with  the  nature  of  the  incoming  flow.  In  most  cases  the  contraction  ratio  is  large  and 
the  flow  is  fairly  slow  but  turbulent.  There  is  generally  no  clear  stagnation  line  at  the  nozzle  intake 
from  which  a  boundary  layer  could  start  laminar.  If  therefore  a  nozzle  throat  is  partly  or  fully  laminar 
this  can  only  result  from  some  form  of  relaminarization  under  a  strong  positive  velocity  gradient. 

In  order  to  discuss  this  problem,  the  criteria  for  relaminarization  due  to  Launder  and  Jones  (ref.26) 
has  been  adopted,  and  applied  to  case  1  (Red  =  3.8  10^)  and  case  2  (Red  =  6.8  10^): 

V  ^  2  5  iQ  5.10^  =>  laminar  flow 
u  dx 


V  ^  <  2.5  to  5.10”^  =>  turbulent  flow 
u  dx 

-  for  throat  Reynolds  numbers  Rcd  =  — - —  well  below  10^,  stays  above  2.5  lO'^  or  even 

jic  u  dx 

5.10'^  from  early  in  the  contraction  down  to  the  close  neighbourhood  of  the  sonic  line.  If,  besides, 
the  surface  is  smooth  there  is  little  doubt  that  the  flow  will  be  laminar  for  a  wide  range  of  upstream 
conditions. 

-  for  high  Reynolds  numbers,  say  well  above  10. 10^,  the  term  does  not  reach  such  high 

u  dx 

values  as  before  and  even  if  the  boundary  layer  is  not  already  turbulent  upstream,  at  the  nozzle  intake, 
it  will  trip  shortly  and  stay  turbulent  till  the  exit. 

-  in  between,  the  situation  is  less  simple  :  assuming  as  above  that  the  boundary  layer  is  initially 

turbulent  (no  stagnation  point !)  the  question  is  whether  will  reach  2.5  10’6  or  5.10‘6  and  stay 

u  dx 

above  such  values  long  enough  and  close  enough  to  the  throat  in  order  to  get  laminar-like  throat 
conditions. 
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Figure  20  collects  the  various  boundary  layer  thickness  parameters  calculated  by  CID  norway, 
whereas  figure  21  depicts  the  value  of  the  Launder  and  Jones  parameter  along  the  nozzle  axis.  As 
figure  21  indicates,  a  relaminarization  of  the  boundary  layer  is  expected  for  toth  cases  and  this  is 
again  visible  in  figure  22. 

Cousteix  and  Aupoix  at  ONERA-CERT  have  observed  that  very  soon  after  ceases  being  high 

u^  dx 

enough,  there  is  a  quick  transition  back  to  turbulent  conditions.  Based  on  such  data,  computations 
have  been  made  for  case  1  (Red  =  3.8  10^)  (fig.23),  showing  again  that  the  boundary  layer  becomes 
lamnar  but  trips  back  to  turbulent  far  enough  from  the  throat  to  get  at  the  sonic  line  a  boundary  layer 
thickness  nearly  equal  to  what  it  would  have  been  if  the  flow  had  never  been  laminar  (fig.24). 

What  may  happen  between  say  0.5  10^  and  about  3.5  10^  is  more  difficult  to  predict  and  may  be 
seriously  depending  on  details  of  the  upstream  conditions.  Especially  whether  it  may  be  considered 
that  some  form  of  stagnation  line  exists  is  certainly  important. 

Basic  experiments  should  be  made  at  these  Reynolds  number  with  various  upstream  geometries  and 
flows,  i.e.  having  or  not  zero  velocity  (stagnation  line)  at  the  upstream  "lips"  of  the  nozzle.  It  will  be 
desirable  to  go  beyond  "single"  back  to  back  comparisons  with  a  reference  nozzle  and  to  include  any 
possible  means  of  boundary  layer  observation.  This  would  help  to  answer  important  practical 
questions  on  the  design  of  the  flowmeters  and  help  to  make  sure  that  flow  rate  measurements  do  add 
properly  when  several  low  Reynolds  nozzles  are  installed  in  parallel. 

A  thorough  understanding  of  these  phenomena  would  also  greatly  simplify  test  operation,  increase 
confidence  and,  as  a  result,  reduce  the  overall  cost  of  these  tedious  calibration  processes. 


33  Real  gas  effects 


Early  in  the  activity  of  the  Working  Group,  it  appeared  that  real  gas  effects  were  not  always 
understood  in  the  same  way  by  the  various  specialists  and  that  tables  and  formulae  used  to  predict 
these  important  effects  should  be  compared. 

In  order  to  avoid  any  misunderstanding  the  equations  allowing  to  perform  such  predictions  and 
comparisons  are  reported  hereafter. 


3.3.1  Notations  (these  notations  apply  to  §§  3.3  through  3.5) 


a  :  Celerity  of  sound  (m/s) 

A  :  Area  (m^) 

B(T)  :  Virial  coefficient 

C(T)  :  Virial  coefficient 

C*  :  Critical  flow  factor  introduced  by  Johnson  (see  equation  (7)) 

C*  ;  critical  flow  factor  for  a  calorically  perfect  gas  =  [Y(2/Y+l)’i'+^/'V“^]^^^ 

ForY=  1.4,  C-  =0.6847314 

Cjjg  :  Correction  factor  for  the  mass  flow  rate  through  a  choked  nozzle  due  to  viscous  effect 

Cdk  :  Correction  factor  for  the  mass  flow  rate  through  a  choked  nozzle  due  to  the  curvature  of 

the  sonic  surface 

:  Correction  factor  for  the  mass  flow  rate  through  a  choked  nozzle  due  to  real  gas  effect 
(virial  effect) :  Cd^  =  Q  real  gas  /  Q  for  a  given  couple  p^, 

Cf^  :  Correction  factor  for  thrust  per  mass  flow  unit  due  to  real  gas  effect 

(virial  effect) :  Ctv  =  |py ^  given  couple  po.  To  -  see  eq.  (18) 

D(T)  :  Virial  coefficient 
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d 

F 

Fv 

h 


Kta 

M 

Mp 

m\ 

M"e 

JC 

p 

Q 

91 

R 


Re 

Re„ 


Re. 


s 

T 

V 

z 


Diameter  of  the  nozzle  throat  (=  2  h)  (m) 

Thrust  (N) 

Thrust  in  vacuum  (or  absolute  thrust)  (N) 

Enthalpy.  Also  radius  of  the  throat  of  a  nozzle  .  \  being  the  geometrical  area  of  the  throat 
:  A,  =  jc  h^ 

Fv 

Absolute  thrust  coefficient :  Kta  =  — - — r 

QfRT^/Q 

Mach  number 

Mach  number  in  the  exit  plane  of  a  contoured  nozzle  (ideal) 

Mach  number  in  the  "inviscid"  part  of  the  flow  at  the  exit  plane  of  a  contoured  nozzle  (see 
eq.  (25  bis)) 

Mach  number  defined  by  eq.  (28  bis).  No  physical  meaning 
Mass  of  one  mole  of  air  (kg) :  =  28.965  x  10'^  kg 

pressure  (N/m^) 

Mass-flow  rate  (kg/s) 

Universal  constant  of  gas  :  =  8.31409  Joule/(mole  .  K) 

21=  8,31409 _ =  287,04  (m/s)2  .  (K)-i  ; 

28,965  X  10-3 

also  radius  of  curvature  of  the  throat  of  a  nozzle  (m) 

Reynolds  number 

Reynolds  number  defined  with  stagnation  conditions  and  half  diameter  of  the  throat : 

r^O 

Reynolds  number  defined  with  critical  conditions  and  diameter  d  of  the  throat : 

j-  Pc  a^  d 

Red  = - 

Pc 

N.B.  :  Re^  =  1.34  Re,,^  for  288  K  stagnation  temperature 
Entropy 

Temperature  in  K  (Kelvin),  °C  (Celsius),  °R  (Rankine) 

N.B.  :  T  (°R)  =  9/5  T  (K) 

Velocity  (m/s) 

Compressibility  factor.  For  a  perfect  gas  :  Z  s  1 


Grecian  notations 


Y 

A 

5(1) 

5(2) 

P 


CO  (M) 

Z(M) 

0(M) 


Ratio  of  specific  heats  (7=  Cp/CJ.  For  a  biatomic  calorically  perfect  gas  :  Y  =  1.4 

Thickness  of  a  boundary  layer 

Displacement  thickness  of  a  boundary  layer 

Momentum  thickness  of  a  boundary  layer 

Viscosity  coefficient  (Poiseuille).  Sutherland's  law  : 

p  =  Po  f  with  T(K) ;  To  =  288  K  ;  C  =  1 10,4  K 

Po  =  1,789  X  10"3  Poiseuille 

Density  (kg/m^) 

-  I  Y-1 
Defined  by  co  (M)  =  1 1  + 


Defined  by  Z  (M)  =  (2  /y+1)2(7-i)  •  ^ ^ 
Defined  by  0  (M)  =  m  (M)  Z  (M)  (l  -h  y 


Y+l 
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Indices 

o  :  Stagnation  conditions 

c  :  Generally  means  critical  conditions  (M  =  1).  But  is  the  geometrical  area  of  the  throat  of 

a  nozzle  (A^  =  k  h^) 

E  ;  In  the  exit  plane  of  a  nozzle 
V  :  In  vacuum 


3.3.2  Real  gas  effects  on  mass  flow  rate 

Frequently  it  happens  that  the  stagnation  pressures  used  in  the  flowmeters  reach  relatively 
high  levels  (10  to  50  bars).  In  these  conditions,  and  at  moderate  stagnation  temperature  levels 
(between  0  to  50°C),  air  does  not  behave  any  longer  as  a  calorically  perfect  gas,  neither  as  a  perfect 
gas.  Then  its  equation  of  state  takes  the  form  : 

p  =  Z  (T,  p)  p  RT 

where  Z  (equal  to  one  for  a  perfect  gas)  is  the  compressibility  factor  known  and  tabulated  by  means  of 
the  virial  coefficients : 


Z  =  1  +  B(T)  p  +  C  (T)  p2  +  D  (T)  p'  +  ...  , 

B(T),  C(T),  D(T)  being  the  virial  coefficients. 

The  entropy-enthalpy  diagram  resulting  from  this  state  equation  allows  to  calculate  isentropic 
expansions  for  various  stagnation  conditions.  More  precisely  : 

-  p„  and  T„  being  a  given  couple  of  stagnation  pressure  and  temperature,  it  is  possible  to  determine 
the  critical  values  (M  =  1)  of  the  density  (p^)  and  of  the  celerity  of  sound  (a^,)  resulting  from  an 
isentropic  expansion.  The  corresponding  mass-flow  rate  for  a  sonic  throat  of  geometrical  area  A  is 
therefore,  for  a  one-dimensional  flow  : 


Q  real  gas  =  p^  .  a,,  .  A^. 


(1) 


-  for  the  same  couple  p„,  T„  of  stagnation  prespre  and  temperature,  it  is  possible  on  the  other  hand 
to  calculate  a  fictitious  mass-flow  rate  assuming  that  air  behaves  as  a  calorically  perfect  gas  (Z  =  1  ; 

Y=  1.4)  all  along  the  isentropic  expansion  from  (p^,  TJ  to  M  =  1. 


This  fictitious  mass-flow  rates  called  Q  ^  for  a  sonic  throat  of  same  geometrical  area  A  as 
previously,  reads,  again  for  a  one  dimensional  flow  : 


Qy=i.4 

Or 

Qy=  1 .4 


Qy=1.4 


Pcy=1.4  •  ^cy=1.4  • 


Ac 


(all  terms  derived  from  po,  Tq 
through  y=1.4  perfect  gas  equations) 


1/2 


Po  Ac 

VRT^ 


21 


Qy  1,4  -  C* 

^  '  VR  lo 

(2) 

(  Y+ 1  ]  1 

(3) 

forY=  1,4,0*  =  0,6847314 

(4) 

It  is  then  possible  to  define  a  coefficient  (V  for  Virial  effect)  which  is  the  ratio  of  the  real 
mass-flow  rate  Q  real  gas  given  by  (1)  to  the  fictitious  mass-flow  rate  Q  y=i,4  given  by  (2) 

(5)  (samep„,T„) 

Qy=i,4 

The  values  of  [C^  -  1]  x  10^  as  a  function  of  p„  and  T„  have  been  calculated  by  Masure  and 

are  given  in  fig.  17  bis-c  extracted  from  [30]  and  in  figure  25.  It  should  be  noticed  that  the  correction 
appears  to  be,  at  constant  T,,,  a  linear  function  of  p„.  Then,  we  can  use,  for  the  domain  of  stagnation 
temperature  and  pressure  considered  here,  the  practical  following  formula  : 


Cdv  =  1  +  0,035 


Po  (atm) 
To(K)  -  210 


(6) 


It  is  worth  noting  that  Masure's  calculations  are  based  on  the  Thermodynamics  tables  of 
Michels,  Wassenaar,  Wolkers  [24]. 


Many  other  authors  made  similar  calculations.  For  instance  Johnson  [11],  combining 
analytical  developments  and  use  of  Thermodynamic  tables  of  Hilsenrath,  Joseph  and  al  [22],  gives  the 
values  of  a  "critical  flow  factor"  C*  defined  as  below  : 


Q  real  g  as 


Pc  •  • 


A,=  C*. 


Po 


Therefore,  (2)  and  (7)  give  : 


(7) 


with  C*  =  0,68473  (see  (4)) 


(8) 


Table  I  gives  some  values  of  as  obtained  by  Masure  (formula  (6))  and  Johnson  (formula 
(8))  for  p  =  10,  20,  30,  40  atm  and  t/=  480°R  (266.7  K)  ;  500°R  (277.8  K)  ;  540°R  (300  K)  ; 
580°R  (322.2  K). 


This  table  leads  to  the  three  following  remarks  : 

a)  The  real  gas  effect  on  mass  flow  rate  is  far  from  negligible.  For  instance  for  p^  =  40  atm  and  T^  = 
266.7  K,  the  correction  due  to  real  gas  effect  is  about  25  x  10'^. 

b)  The  actual  mass-flow  rate  is  always  higher  than  the  mass-flow  rate  calculated,  were  the  air  a 

calorically  perfect  gas  (Z  =  1  ;y  =  1.4). 

c)  Values  given  by  Masure  and  Johnson  are  very  near  each  other  but  different  :  the  difference 
between  them  never  exceeds  0.6  x  10'^  on  table  I.  The  ratio  of  these  two  values  is  presented  on 
the  figure  26  for  a  larger  range  of  pressure  :  0  <  p^  <  60  bar.  On  this  figure,  the  gap  is  slightly 
smaller  than  1.10'^.  Since  Masure  and  Johnson  made  what  is  indeed  the  same  calculation  in  two 
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different  ways  but  starting  from  two  different  sources  of  thermodynamic  data,  reducing  the  gap 
would  need  to  recheck  thermodynamic  tables  ! 


Remark :  the  pressure  range  for  formula  (6)  is  :  0  atm  to  about  50  atm.  For  higher  pressures 
(50  ^  Po  <  100  atm),  it  is  recommended  to  use  Johnson's  values  of  the  critical  flow 
factor  C*  given  in  [1 1]  and  to  apply  equation  (8). 

Conclusions 


If  high  accuracy  is  needed  when  measuring  mass  flow  rate  with  sonic  Venturis,  real  gas 
effect  (virial  effect)  must  be  taken  into  account. 

The  actual  mass  flow  rate  is : 


q  real  gas  =  .  C*  .  (9) 

Viv  Iq 

Were  air  a  calorically  perfect  gas  with  Y=  1.4,  then  is  equal  to  one.  The  difference  - 
1  can  be  substantial. 

Refining  real  gas  corrections  (Cdv)  beyond  1.10'^  would  require  to  check  thermodynamic 

tables. 


Table  I 

Real  gas  effect  (virial  effect)  on  mass  flow  rate  for  a  sonic  Venturi 


T 

(atm) 

O 

10 

20 

30 

40 

(°R) 

t„(kK 

(J) 

1  +  6.1  X 

10’^ 

(J) 

1  +  11.9  X  10'^ 

(J) 

1  +  17.9  X 

10' 

(J) 

1  +  24.1  .  10' 

480 

266.7 

(M) 

1  +  6.2  X 

10-5 

(M) 

1  +  12.3  X  10'^ 

(M) 

1  +  18.5  X 

10' 

(M) 

1  +  24.7  X  10  ' 

(J) 

1  +  5.4  X 

10-^ 

(J) 

1  +  10.5  X  10'^ 

(J) 

1  +  15.6  X 

10' 

(J) 

1  +  20.8  X  10  ' 

500 

277.8 

(M) 

1  +  5.2  X 

10-' 

(M) 

1  +  10.3  X  10-^ 

(M) 

1  +  15.5  . 

10' 

(M) 

1  +  20.7  X  10  ' 

(J) 

1  +  4.0  X 

10^ 

(J) 

1  +  8.0  X  10  ^ 

(J) 

1  +  11.9  X 

10' 

(J) 

1  +  15.7  X  10  ' 

540 

300 

(M) 

1  +  3.9  X 

10-^ 

(M) 

1  +  7.8  X  10-3 

(M) 

1  +  11.7  X 

10' 

(M) 

1  +  15.6  X  10-' 

(J) 

1  +  3.0  X 

10-" 

(J) 

1  +  6.1  X  lO'^ 

(J) 

1  +  9.0  X 

10' 

(J) 

1  +  11.9  X  10  ' 

580 

322.2 

(M) 

1  +  3.1  X 

lo-^ 

(M) 

1  +  6.2  X  10-5 

(M) 

1  +  9.4  X 

10' 

(M) 

1  +  12.5  X  10  ' 

Coefficient  C^^  =  ^^fll^for  a  sonic  venturi  as  given  by  Johnson  (J)  and 


Masure  (M) 


Remark  :  Johnson's  C^^  values  have  been  deduced  from  data  given  in  table  II  of  reference  [1 1]  using 
equation  (8)  of  the  present  report. 

Masure's  Cj^  values  have  been  deduced  from  equation  (6)  of  the  present  report. 
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3.3.3  Real  gas  effects  on  thrust 


For  a  subsonic- supersonic  nozzle,  the 
thrust  in  vacuum  Fy  (also  called 
absolute  thrust)  can  be  calculated  using 
momentum  equation  : 

Fv  =[(pe  +Pe  Ve)Ae]  X 


where  x  is  a  unitary  horizontal  vector 
oriented  upstream. 


The  quantity  (pE  +Pe  Vg)  Ae  x  is  equal  to  the  integral  of  pressure  on  the  internal  surface  a  b  c 
d  e  f  g  h  i  j  k  of  the  nozzle. 

We  now  calculate  F.^  assuming  one-dimensional  isentropic  flow  for  both  cases  : 

1)  air  is  assumed  to  be  a  perfect  gas  (compressibility  factor  equal  to  one  :  Z  =  1)  with 
constant  specific  heats  (y=  1-4). 


2)  air  is  considered  as  a  real  gas  (Z^  1). 


1  Case  1  :  air  is  assumed  to  be  a  caloricallv  perfect  gas  (Z  =  1  :  y  =  1 .4) 

Fv=  ',Pe Pe  ^'eI  Ac 


—  PeU  -i- yVIc)  Ac  (10) 

(E  means  :  in  the  exit  plane  Ag,  Mg  :  Mach  number  in  Ag,  A^  :  sonic  throat  area,  P„  :  stagnation 
pressure  of  the  upstream  flow). 

On  the  other  side,  in  the  frame  of  the  same  assumption,  the  mass-flow  rate  is  : 


Q 


wi'h  C*  =  } 

fory=  1.4  :C*  =  0,6847314 
Tq  ;  stagnation  temperature  of 
the  upstream  flow  (K) 


Combining  (10)  and  (11)  gives  : 


Or 


QVin^/c*  ^ 


=  Ei.:^.(l+yM|) 


(11) 


(12) 
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In  order  to  recall  that  the  ratio  q  /  C*  ^Qual  to  the  right-hand  side  of  (12)  only  if  air 
is  assumed  to  be  a  calorically  perfect  gas,  we  use  hereafter  the  abbreviation  y  =  1 .4  : 


(J„.,,4,yRT7/Ci  P»  A: 


(13) 


Remarks : 


Pe. 

Po 

Al 

Ac 


/  Y-  1  ..'r'l'/'''-! 

=  I  +  ^Mif 


!  -) 


1 


y-1) 


R  =  defined  in  3.3.1  as  the  ratio  of  the  universal  constant  of  gases  to  the  molar  mass  of  air. 


2  Case  2  :  air  is  considered  as  real  gas  (Z^  1) 


being  now  the  thrust  in  vacuum  in  real  gas  conditions,  we  have  : 


Fv=  (pe  +  PeVe)Ae 

=  Pe  Ae  +  (Pe  Ve  Ae)  Ve 


(all  quantities  consistent  with 
isentropic  real  gas  expansion) 


(14) 


Starting  from  given  values  of  the  stagnation  pressure  p„  and  stagnation  temperature  T„,  it  is 
possible,  using  either  thermodynamic  tables  of  air  or  alytical  formulas  taking  into  account  the  virial 
effect  for  air,  to  calculate  isentropic  expansions  of  the  gas.  For  a  particular  value  of  the  pressure, 
called  critical  pressure  p^,  the  velocity  of  the  flow  is  equal  to  the  celerity  of  sound  :  V  =  a,.  =  and  the 
corresponding  value  of  the  density  is  p^  If  A^  is  the  area  of  the  throat  (assumed  sonic),  the  mass-flow 
rate  is  ; 


Q  real  air  =  .  a^ .  A^, 

For  any  other  pressure  p^  below  p^,  the  values  of  the  velocity  of  the  flow,  of  the  celerity  of 
sound  and  of  the  density,  for  the  same  isentropic  expansion,  will  be  called  Vg,  a^,  pg. 

To  calculate  the  area  Ag  occupied  by  the  flow  in  this  new  state,  we  use  the  continuity 
equation : 

Pe  •  ^E  •  Ae  =  Pc  a^  Ap  =  Q  real  air 


This  gives  : 


Q  real  air 

Pe  •  ^E 


(15) 


Therefore,  (14)  becomes  : 


Vreal  air 


^  Q  real  air  ^  .  ,, 

=  Pe-  ^ - ^+9realair.  Vg 

Pe-  Ve 


(16) 


We  now  introduce,  in  the  same  way  as  for  case  1,  the  ratio  : 


Vreal  air 


Q  real  air  YR  l  p  /  C* 


Using  (16)  we  find  : 

_Ze_+ Ve 

FVrcal.ir  _PeVe 

Qreal  air  ^  O 


(17) 


It  is  now  interesting  to  compare  the  two  ratios  as  given  by  equations  (13)  and  (17)  for  the 

same  nozzle  defined  by  A^/A^,  each  of  these  two  ratios  being  in  fact,  besides  the  coefficient  VRTq/Q  ^ 
the  unitary  thrust  (that  is  to  say  the  thrust  per  mass  flow  unit)  for  each  of  the  two  cases  1  and  2  . 


We  intt-oduce  consequently  a  coefficient  Cj.^  (T  for  thrust,  V  for  virial)  defined  by  : 


^  _  Fvreal  ai/Qrcal  air 

^  Fv(^l,4/Q(')fcl,4) 


(18) 


Masure  [30],  carrying  out  calculations  using  thermodynamic  tables  of  reference  [24]  in  the 

domain  0  <  p,,  <  40  atm,  0°C  <  To  <  75°C  and  for  different  values  of  Ag/A^  (Ag/A^  =  1  ;  1.7  ;  3  ;  oo 
(ultimate  expansion)),  found  that,  for  a  given  nozzle  and  a  given  stagnation  temperature  To,  the 
coefficient  is  a  function  practically  linear  of  the  stagnation  pressure  p„  (like  for  the  mass  flow 
rates).  The  effect  of  pressure  may  therefore  be  outlined  by  putting  it  in  the  form  indicated  in  fig.  17 
bis-d,  where  the  results  of  the  calculations  are  presented  for  various  A^/A^,. 

Let  us  recall  that  the  statements  made  here  for  the  thrust  are  only  valid  within  a  domain  of 
stagnation  pressure  identical  to  that  considered  for  mass  flow  rates. 

Let  us  consider  an  example  :  p„  =  40  atm ;  T^  =  0°C  ;  Ag/A^.  =  3 

=  1  -  24.4  X  10-3 

Fv(^,,4/Q(T=1,4) 

Such  a  correction  is  far  from  negligible. 

Johnson  [11]  carried  out  similar  calculations  which  can  be  compared  to  Masure’s  as  it  has 
been  done  for  thrust  in  above  3.3.2. 

Such  a  comparison  will  be  carried  out  for  sonic  nozzle  only  because,  Johnson  did  not  carry 
out  his  isentropic  expansions  beyond  critical  flow  (M  =  1) 

For  a  sonic  nozzle,  equation  (17)  reads  : 

F  P  ^ 

Q,eai»vinyc*  VRVc* 
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equation  (19)  becomes : 


^real  air 


Pc/Po 


+  G 


QxealairVRT^/C*  VRT^ 


(20) 


Johnson's  results  [11]  can  now  be  used  to  calculate  the  right-hand  side  of  (20).  Choosing 
values  for  the  stagnation  pressure  and  stagnation  temperature  T„,  one  finds  in  [1 1]: 

a)  the  value  of  p^p^  (see  table  11(c) :  Critical  pressure  ratio) 


(3)  the  value  of  factor  table  jj  (a) :  critical  flow  factor) 


y)  the  value  of  a, :  a,  is  equal  to  the  nozzle-throat  velocity  because  M  =  1.  In  table  11(b), 
Johnson  gives  the  values  of  the  ratio  of  nozzle-throat  velocity  to  the 
speed  of  sound  at  1  atmosphere  and  491.688°R  (1087.42  ft/sec).  Remark 
:  1087.42  ft/sec  =  331.4456  m/s. 

With  these  values  and  with  given  bv  f4)  and  R  given  following  3.3. 1  by 

R  =  8.31409, ^  287.04  (m/s)^  (K)-^ 

28.965x10-" 

Fv 

the  ratio  - — -  can  be  determined  thanks  to  (20). 

Qreal  air  V  RTo/Cj 

The  value  of  the  same  ratio  for  7=  1.4  is  given  by  (13)  in  which  one  has  to  assume  that  the 
nozzle  is  sonic  (M^  =  1) : 


Fvi„.i 


Qri.4  m7/cl 


1.26788 


Finally,  we  obtain  Ct.  = 

Fvy=iVQi-i.4  1.26788 


(21) 


Johnson's  Cj.  values  as  given  by  (21)  have  been  calculated  for  p^  =  10  ;  20  ;  30  ;  40  ;  50 
atm  and  To  =  480°R  (266.7  K) ;  500°R  (277  K) ;  540°R  (300K) ;  580°R  (322.2K) ;  620°R  (344.4K). 

Results  are  presented  on  figure  27.  It  appears  clearly  that  Cj.  is  almost  a  linear  function  of  p^ 
for  a  given  T^,  confirming  Masure's  conclusion. 

The  comparison  between  Johnson's  results  and  Masure's  results  concerning  sonic  nnrvJp.  is' 
presented  on  table  11  and  figure  28. 

Figures  27  and  28,  and  table  II  lead  to  the  following  remarks  : 
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a)  The  real  gas  effect  on  the  coefficient  Cr  is  far  from  negligible  for  a  sonic  nozzle.  For 
instance  for  =  50  atm  and  T„  =  266.7  K,  the  correction  due  to  real  gas  effect  is  about  - 
24  X  lO  l 

b)  The  actual  unitary  thrust  (thrust  per  unit  of  mass  flow  rate)  is  always  lower  than  the 
unitary  thrust  calculated,  were  air  a  calorically  perfect  gas  (Z  =  1 ;  y  =  1-4). 

c)  Values  given  by  Masure  and  Johnson  are  very  near  each  other  :  the  difference  between 
them  seldom  exceeds  1.10'^  (see  table  II  and  figure  28). 

For  a  supersonic  nozzle,  figure  17  bis.d  shows  that  :  the  higher  Ag/A^.,  the  higher  the  real  gas  effect 
on  thrust. 


Table  II 

Real  gas  effect  (virial  effect)  on  thrust  for  a  sonic  nozzle 


T. 

(°R) 

\,^(atm) 

t„(kN. 

10 

20 

30 

40 

50 

(J) 

1  -  4.9  X  10’^ 

(J) 

1  -  9.8  X  10'^ 

(J) 

1  -  14.4  X  10’’ 

(J)  1  -  19.1  X  10’ 

(J)  1 

-  23.4  X  10  ’ 

480 

266.7 

(M) 

1  -  4.7  X  10-^ 

(M) 

1  -  9.5  X  lO'^ 

(M) 

1  -  14.2  X  10  ’ 

(M)  1  -  18.9  X  10  ’ 

(M)  1 

-  23.6  X  10  ’ 

(J) 

1  -  4.3  X  10-^ 

(J) 

1  -  8.3  X  10-^ 

(J) 

1  -  12.2  X  10  ’ 

(J)  1  -  16.0  X  10  ’ 

(J)  1 

-  19.7  X  10  ’ 

500 

277.8 

(M) 

1  -  4.2  X  10-' 

(M) 

1  -  8.4  X  10’^ 

(M) 

1  -  12.6  X  10  ’ 

(M)  1  -  16.8  X  10  ’ 

(M)  1 

-  21.0  X  10’ 

(J) 

1  -  3.0  X  lO'" 

(J) 

1  -  5.9  X  10  ’ 

(J) 

1  -  8.7  X  10  ’ 

(J)  1-11.4x10’ 

(J)  1 

-  13.8  X  10  ’ 

540 

300 

(M) 

1  -  3.0  X  10-^ 

(M) 

1  -  6.1  X  10  ’ 

(M) 

1  -9.1  X  10’ 

(M)  1  -  12.2  X  10  ’ 

(M)  1 

-  15.2  X  10  ’ 

(J) 

1  -  2.1  X  10-^ 

(J) 

1  -4.1  X  10’ 

(J) 

1  -  6.0  X  10  ’ 

(J)  1  -  7.8  X  10  ’ 

(J) 

1  -  9.6  X  10  ’ 

580 

322.2 

(M) 

1  -  2.1  X  10^ 

(M) 

1  -  4.2  X  10  ’ 

(M) 

1  -  6.2  X  10  ’ 

(M)  1  -  8.3  X  10  ’ 

(M) 

-  10.4  X  10  ’ 

(J) 

1  -  1.3  X  lO'" 

(J) 

1  -  2.7  X  10  ’ 

(J) 

1  -  4.0  X  10  ’ 

(J)  1  -  5.2  X  10-’ 

(J) 

1  -  6.2  X  10  ’ 

620 

344.4 

(M) 

1  -  1.3  X  10-" 

(M) 

1  -  2.5  X  10  ’ 

(M) 

1  -  3.8  X  10  ’ 

(M)  1  -  5.1  X  10  ’ 

(M) 

1  -  6.3  X  10  ’ 

Coefficient  Cpy  =  fQj.  ^  nozzle  as  given  by  Johnson  (J)  and  Masure  (M) 


Remark :  Johnson's  Cj  values  have  been  deduced  from  data  given  in  table  II  of  reference  [11] 
using  equations  (20)  and  (21)  of  the  present  report. 

Masure's  Cj.^  values  are  deduced  from  figure  17  bis-d. 


28 


3.4  Viscous  effects  on  thrust 


3.4.1  Preliminary  remarks 

Addressing  the  general  issue  of  viscous  effects  on  thrust  would  obviously  be  far  more 
complicated  than  doing  it  as  above  for  mass  flow  rates.  Indeed,  amongst  others,  the  following 
difficulties  would  be  met : 

•  boundary  layer  effects  of  all  kinds, 

•  shock  boundary  layer  interaction 

•  computation  of  the  supersonic  free  flow  in  presence  of  a  boundary  layer 

•  separation  etc... 

Nevertheless,  it  seemed  useful,  within  the  scope  of  this  report  to  present  some  relatively 
simple  calculations  which  provide  in  most  practical  cases  a  reasonably  good  approximation  of  thrust 
coefficients.  This  should  allow  operators  of  thrust  measuring  rigs  to  do  some  first  cross-checking 
with  their  force  measurements.  The  simplifying  physical  assumptions  are  the  following  : 

•  thin  boundary  layers  at  the  throat  and  in  the  exit  plane 

•  no  separation  in  the  nozzle, 

•  no  dramatic  effect  of  the  boundary  layer  in  the  nozzles  on  the  isentropic  expansion  of  the 
supersonic  free  flow  :  the  boundary  layer  changes  slightly  the  «  shape  »  of  the  contoured  nozzle, 
hence  changing  the  exit  Mach  number  but  the  free  flow  at  the  exit  plane  remains  reasonably  uniform. 

Calculation  presented  hereunder  are  made  under  these  assumptions,  following  [31] : 

3.4.2  Non  viscous  thrust  reference 


To  estimate  viscous  effects  on  thrust  we  assume  in  this  §  3.4  that  air  is  a  calorically  perfect 
gas  (Z  =  1  ;  Y  =  1.4). 


Let  us  consider  an  axisymmetric  contoured  subsonic-supersonic  nozzle.  The  shape  of  this 
nozzle  has  been  calculated  in  such  a  way  that  the  flow  in  the  exit  plane  (E)  of  the  nozzle  is  supersonic, 
uniform  and  parallel  to  the  axis  x'x  of  the  nozzle,  the  viscous  effects  being  ignored  for  the  moment.  In 
the  vicinity  of  the  throat,  the  contour  is  assumed  to  be  circular  (radius  of  curvature  R).  The  ratio  R/h, 
h  being  the  radius  of  the  throat,  is  assumed  to  be  in  the  order  of  what  it  uses  to  be  on  typical  such 
nozzles  i.e.  about  4. 

being  the  thrust  in  vacuum,  \  the  area  of  the  geometrical  throat  (A^  =  Tih^),  p,,  p,,  M.  the 
density,  pressure  and  Mach  number  in  (E)  respectively,  p^  and  T  the  stagnation  pressure  and 
temperature,  we  have  (see  (10)) : 
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Fv  =  ^  .  (i+yMe)  •  Po  Ac  with  Ae  =  tiRe  (10) 

Po  PiQ 

Due  to  non  uniformity  of  the  flow  in  the  plane  of  the  geometrical  throat,  the  area  of  the 
(fictitious)  critical  throat  is  not  but  where  0)^^  is  the  coefficient  taking  into  account  the 

curvature  of  the  sonic  surface  <  1).  The  mass  flow  rate  is  therefore  : 


Q  =  Cd, 


one  dimensional 
flow 


)  =  Cdk  •  C 


* 

i 


Po  Ac 


We  define  now  an  absolute  (i.e.  in  vacuum)  thrust  coefficient  as  : 


Kta  = 


Fy 

QVRTo/C- 


(22) 


(23) 


(22)  shows  that : 


and  (10)  shows  that : 


Kta  = 


Fy 

Po  Cdk  Ac 


Kta  = 


PE 
Po  ■ 


Ae 

Ac  Cdk 


(i+ymI) 


(24) 


_  _  /  Y_j  Vy/y-i 

Now  :  pe/Po  =  (o(Me)  with  co(M)  =  [l  + 

Ab/AcCdk  =  2:(Me)  (24bis)  with  E(M)  =  [l  + 

so,  Kta  =  (KMe)  with  (t)(M))  =  co(M)E(M)(l +yM^)  (25) 

(inyiscid  flow  ;  y  =  cte) 

For  instance,  for  Y  =  1.4  and  Me  =  2,  K-pA  =  1,42342 
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3.4.3  A  simplified  method  to  estimate  viscous  effects 


Along  the  wall  of  the  nozzle  a  boundary  layer  appears,  the  effect  of  which  is  to  reduce  the 
mass  flow  rate,  for  given  p^.  To,  and  to  modify  the  mean  Mach  number  in  the  exit  plane  (E).  The 
thrust  coefficient  defined  by  (23)  will  therefore  be  somewhat  different  from  that  given  by  (25). 


Ag  =  7c  Rg 


(E)^ 

Fv  being  the  new  thrust  in  vacuum,  the  momentum  equation  permits  to  write  : 


(P 

0 


Fv  =  (p+pV^l  27rydy 


Let  Mg'  be  defined  by  e(M£)  =  — ^ttRe  Se 

Ac  Cdk  •  Cds 


(25  bis) 


where  is  the  coefficient  of  mass  flow  rate  reduction  due  to  the  boundary  layer  at  the 
x(i) 

throat  of  the  nozzle  and  where  Oe  is  the  displacement  thickness  of  the  boundary  layer  in  the  exit  plane 
(E). 


Fv  =  Pe  Ae  +  I  pV^  .  2  Tiydy  ,  the  pressure  p'g  ,  assumed  uniform  in  the  exit  plane  (E), 

Jo 

being  connected  to  the  upstream  stagnation  pressure  p„  through  : 

Pe/Po  =  ^Mg)  |we  recall  that  co(M)  =  1 1  +  ^  M^j  ^  ^  | 

Let  pg  and  Vg  be  the  values  of  density  and  velocity  of  the  flow  in  the  exit  plane  (E),  outside 
of  the  boundary  layer  :  these  values  are  connected  to  the  upstream  stagnation  conditions  through  the 


above  mentioned  Mach  number  M'g. 


PqAc  )  . 


The  new  mass  flow  rate  Q',  now  equal  to  Cds  |Cdk  •  Cj  |  is  also  given  by  : 

Q'  =  PEVEJc(RE-5l'f  < 


Moreover,  taking  into  account  the  definition  of  the  momentum  thickness  Se  ^  in  the  exit  plane 
(E),  that  is  : 
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[  pv(ve-v) 

Jr-a 


27iydy  =  PeVe  27iRe  5e  •  Ve  , 


it  is  easy  to  demonstrate  the  following  relation  : 

Fv  =  Pe  .  Ae  +  Q'  Ve  -  2k  Re  5®  Pe  Ve^  (27) 

Let  us  show  now  that,  following  [32],  the  relation  (27)  can  be  simplified.  Consider  a 
fictitious  nozzle  fed  with  a  non  viscous  gas  of  same  stagnation  pressure  and  temperature  p,,,  T„,  as 
the  real  one,  having  the  same  exit  area  Ag  and  the  same  mass  flow  rate  Q'  as  the  real  one. 

Since,  by  hypothesis,  p„,  T„,  Q'  are  the  same  in  the  two  nozzles,  the  critical  area  A,^  of  the 
fictitious  nozzle  is  given  by  : 

Ks  =  K-  Cd^  •  Ccg 

Let  Fv  be  the  thrust  in  vacuum  of  this  fictitious  nozzle  fed  with  air  considered  as  a  non 


viscous  gas : 


Fv  -  Pe  Ae  +  Q'  Ve 
and  = 


(28  bis) 


From  (27)  and  (28) : 


Fv  -  Fv  =  PeAe  +  Q  Ve  -  2jrRE  5e  Ve^e  '  (pe-^e  +  Q  Ve)  (29) 

As  Q'  =  Pe  Ve  Ae,  (29)  reads : 

F;  -  F;  =  Ae[(pe  -  Pe)  +  Pe  V^Ve  -  V^]  -  2jtRE  Sg’  pj  V'i  (30) 

But  p'g  and  p''^  can  be  considered  as  the  values  of  the  pressure  in  two  neighbouring  states 
during  an  isentropic  expansion  from  the  stagnation  pressure  p„,  these  two  neighbouring  states  being 
defined  by  the  Mach  numbers  M'^  and  M"g,  near  each  other.  It  is  exactly  the  same  for  the  two 
velocities  V'^  and  V  g. 

I  Pe  ■  Pe  ”  ^Pe 

We  may  therefore  write  :  { 

lVE-VE  =  dVE 

On  the  other  hand,  h  being  the  enthalpy  and  s  the  entropy,  we  have  ; 


and,  as  ds  =  0, 


dh  =  Tds  +  dp/p 
dh  +  VdV  =  0 

p  VdV  +  dp  =  0 


(thermodynamics) 

(first  law  of  thermodynamics) 


Consequently  :  Pe  '  Pe  +  Pe  ^e  (Ve  -  Ve)  =  0  and  (30)  reads : 

f;  -  f;  =  -  27tRE  .  5^  .  Pe  .  w'i 


32 


Or,  with  an  error  of  second  order : 


f;  =  -  2  tiRe  .  ^ .  p'e  .  (31) 

Dividing  the  left-hand  side  of  (31)  by  Q  YRTo/C*  and  the  right-hand  side  of  the  same 
equation  by  ,  one  obtains : 


0{M"e)  -  TCo(M"e)  X(M"e)  M’i  .  ^ 

Re 


{actual  nozzle 
viscous  flow 
7=1.4 


We  recall  that  (t)(M)  =  co  s(l  +  yM^)  (see  (25))  and  that  the  Mach  number  M  ^  is  the 
(supersonic)  solution  of  x(Me)  =  Ae^Ac-Cdk-Cds)-  This  Mach  number  has  no  physical  meaning.  One 

s;(i) 

sees  that,  to  apply  (32),  it  is  necessary  to  know  the  displacement  thickness  be  at  the  throat  of  the 
nozzle  because  Cdj  =1  — and  the  momentum  thickness  Sg  ^  in  the  exit  plane  (E). 

Use  of  (32)  will  be  presented  later  (see  §  3.5). 

Another  writing  of  equation  (321 

It  is  possible  to  give  another  formulation  of  equation  (32)  in  which  not  only  5^  ^  but 
also  5c^  appear  explicitly. 

We  first  remember  that  the  Mach  numbers  Mg  and  Mg,  which  are  near  each  other,  are 
the  supersonic  solutions  of  the  following  equations  : 


S  (Me)=x^ 

^  ^C-'-Dk 


(see  24(bis)) 


Ac.Ce^.Cds 


(see  28(bis)) 


Then,  neglecting  second  order  terms  : 


(1)(Me)  =  (KMe)  +-%-{Me)  [Z  (Me)  -  S  (Me)] 

dZ. 

As  4)(M),  ^  (M),  co(M)  verify  =  CO,  it  is  easy  to  obtain  after  some  straight¬ 


forward  but  lengthy  calculations  : 


viscous  flow 
7=  1,4 


<KMe)  |l  +  .  2  wj 

1  l+yM|  l-fyMi  } 
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Or: 


Kta, 


viscous  flow  I 
7  =  1.4  i 


=  Kta 


non  viscous  flow! 
7  =  1.4  / 


.1^ 


25^^/h 


l+yMl  l+yM^ 


YM| 


2  5^2^) 

l/ 


(32  bis) 


35  Practical  formulae  to  calculate  mass  flow  rate  and  thrust 

3.5.1  Introduction 

In  §§  3.3  and  3.4  several  formulae  have  been  proposed  to  calculate  mass  flow  rates  and 
thrust  coefficients.  It  is  proposed  now  to  show  how  such  formulae  can  be  used  practically.  This  will 
allow  to  give  examples  of  the  relative  magnitude  of  the  various  correction  factors. 

3.5.2  Practical  formulae  to  calculate  mass  flow  rate  in  a  sonic  Venturi 


The  actual  mass  flow  rate  is  : 

Q  =  Cd..Cd..Cd.C*;^  (33) 

with  : 

p„ :  stagnation  pressure  of  the  upstream  flow  (N/m^) 
To  :  stagnation  temperature  (K) 

R  =  287.04  (m/s)'  (K)'" 


Ac  :  Area  of  the  geometrical  throat  (Tch')  (m') 
assumed  choked 


C*  =0.6847314 

Cn  :  coefficient  due  to  the  curvature  of  the  sonic  surface 

Cd  :  coefficient  due  to  boundary  layer  effect 

1  ^(i)\ 

(CDa=l-2^) 

:  coefficient  due  to  real  gas  effect. 

As  an  example,  let  us  assume  that  the  sonic  Venturi  is  the  ONERA  short  radius  nozzle  R^d^ 
=  2  (see  appendix  F)  with  d^,  =  20  mm,  p^  =  45  bar.  To  =  288  K.  We  have  now  to  determine  Cq  , 

K 

c  c 

'-'dg’  '-Dy- 


a)  For  d^/R^  =  0.5,  figure  17  bis-a  gives  : 
Cj.  =0.99857 


(34) 
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P)  Cog  For  dyR^  =  0.5,  assuming  a  turbulent  boundary  layer  at  the  throat  and  following  simple 
boundary  layer  calculations  used  in  ref  [23]  and  in  3.2.2.3,  figure  17bis-b  gives  : 

{l-CD.)Ret(5=  0.0454  (35) 

with  Reo.h  =  ^  ^  =  10,45  x  10^ 

M^o 

=>  Cd  =0.99693 

Y)  Cdv  For  p„  =  45  bar  (that  is  for  p„  =  44.41  atm)  and  =  288  K,  figure  17bis-c  gives  the  value  of 
(Cj)^  -  1)  X  10^  after  interpolation.  It  is  easier  to  use  the  formula  : 

^°v=> -0.035  (see  (6)) 

we  obtain : 

Cdv  =  1  +  19. 

Summary : 


The  actual  mass  flow  rate  is  therefore,  from  (33) : 


93  X  10-3  =  1.01993 

Cd  =  0.99857  =  1  -  1.43  X  lO'" 

‘Ok 

Cd  =  0.99693  =  1  -  3.07  x  10^ 

=  1.01993  =  1  +  19.93  x  10’^ 


Q  =  (1.01993)  ( 0.99693)  (0.99857)  (0.6847314)  45  x 


Jt(2.10-f/4 
■  ^287.04  X  288 


3.4184  kg/s 


3.5.3  Practical  formulae  to  calculate  thrust  coefficient  Kj^. 
Another  nozzle  will  be  considered  for  this  example. 

F 

Let  us  recall  that  Kr,  =  — ; — ^ — - 

Q-Woic] 


We  will  calculate  for  the  ONERA  =  2  supersonic  nozzle.  Let  us  recall  that  The 
throat  diameter  is 

d^  =  61.64  mm  and  the  diameter  in  the  exit  plane  (E)  is  :  2  =  80,0154  mm. 

The  shape  of  this  axisymetric  nozzle  has  been  calculated  as  it  is  explained  at  the  beginning  of 
§  3.4.  In  particular  the  radius  of  curvature  of  the  nozzle  in  the  vicinity  of  the  throat  is  R  =  123.28  mm, 
that  is  2  times  the  throat  diameter. 
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a)  Hypothesis  n°  1 :  Inviscid  flow  and  air  is  a  calorically  perfect  gas  (Z=  1  ;  y  =  1,4) 
Equation  (25)  gives  :  =  0(Me)  with  =  2 

Hence ; 


Kta 

=  1.42342 

(inviscid  flow) 

(7  =  1.4 

) 

(37) 


P)  Hypothesis  n°  2 :  Viscous  flow  (turbulent  boundary  layer)  and  air  is  a  calorically  perfect  gas  (Z  = 
l;y=1.4) 


To  apply  (32  bis),  we  must  calculate  5^^’^  and  (or  measure). 

Let  us  assume  p^  =  8  atm  ;  =  288  K 

|'Reo,h  =  5.7  X  106 
Then  | 

1(1  -  Cd,)  Re;;6  =  0. 0454  (fig  17  bis-b)  (38) 


From  (38) :  Cq  =  0.9966  and  ^/h  =  1.7  x  lO'^ 

In  the  exit  plane,  the  value  of  can  be  obtained  through  two  ways  :  computation  (e.g.) 
integration  of  the  integral  equation  of  Von  Karman)  or  direct  measurement  in  the  exit  plane  with  a 
Pitot  rake. 


According  to  [33],  the  results  are  : 

Integration  of  integral  equation  of  Von  Karman 


s-(l) 

5e  _  14.2 
h  1000 


er 


4.74 

1000 


=  3 


Direct  measurement : 


Let  us  retain  the  first  value  for 


18.8 

1000 


5.1 

1000 


5^ 


=  3.7 


5e  ^  ^  4.74  _  5e  ^  h  ^  4.74  61.64  ^  3,65 

h  1000  Re  h  Re  1000  80.01  1000 


Then,  using  (32  bis) : 

=  1.41534 

(viscous  flow] 

)y=L4  I 
ipo  =  8atm 
(To  =  288  K  ) 


(39) 
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Let  us  assume  now  that,  the  temperature  T„  being  unchanged  (T^  =  288K),  the  pressure  p„ 
is:  p„  =  48  atm.  (i.e.  six  times  the  former  pressure) 

Re,,^  =  (5.7  X  10")  X  6  =  34.2  x  10" 

(38)  =i>  (l -  =  0.0454  ^  =  1  -  2.50  x  10“'  =  0.9975  //z  =  1.25  x  10"' 

_  2.71 

Re  ~  1000 


Then,  using  (32  bis) : 


Y)  Hypothesis  n°  3  :  Viscous  flow  (turbulent  boundary  layer)  and  air  is  no  longer  a  perfect  gas  (Z^l) 
For  A/A^  =  1.7  and  To  =  288  K  (s  15°C),  fig  17  bis.d  gives  : 

I^TAreal  air  _  J  ^  X  10^  —  _  4  35 

lKTAr=i.4  f  Po(atm) 

If  Po  =  8  bar  (=  7.90  atm),  ^TAr^ai  air  =  j  _  3  44  x  10-3  =  0.99656 

^TAr=i.4 

If  Po  =  48  bar  (=  47.37  atm),  1  -  20.61  x  10-3  =  0.97939 

^TA=1.4 

Using  the  values  given  by  (39)  and  (40),  one  obtains  : 
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Summary  :  Values  of  obtained  in  a ,  p  ,  y  are  summarized  beneath 


a  inviscid  flow 
perfect  gas 
(Z=l;y=1.4) 

-> 

1.42342 

n  viscous  flow  j 

'^perfect  gas  j 

(-Z=l;y=1.4l}  t 

fPo  =  8  atm 

-> 

1.41534 

To  =  288  K  1 

[po  =  48  atm 

-> 

1.41741 

y  viscous  flow 
real  gas 

(z^  l]  * 

|Po=  8  atm 

1 

-> 

1.41047 

To  =  288  K 

^o  =  48  atm 

-> 

1.38820 

This  example  shows  that,  when  measuring  the  thrust  coefficient  of  a  nozzle  with  a  bench,  attempts  to 
cross-check  experimental  data  with  theoretical  predictions  must  take  into  account  both  viscous  effects 
and  real  gas  effects  in  order  to  have  any  meaning. 
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4.  EXAMPLES  OF  EXPERIMENTAL  RESULTS 
4.1  Boeing's  experience  on  ASME  and  "cubic"  nozzles 


The  ASME  nozzle  (fig.  14)  has  long  been  used  as  a  standard  thrust  nozzle  in  many  wind  tunnel  and 
propulsion  laboratory  facilities.  As  accuracy  of  mass  flow  and  thrust  measurements  has  improved, 
discrepancies  have  been  noted  in  the  performance  of  the  ASME  nozzles  when  compared  to  long 
accepted  characteristics.  Reference  10  was  one  of  the  early  publications  to  define  the  source  of  the 
ASlSffi  nozzle  problems. 

A  new  nozzle,  called  "cubic  nozzle"  was  designed  in  1978  by  Boeing  Propulsion  Laboratory  to  have 
improved  characteristics  as  the  result  of  a  thinner  boundary  layer  and  no  cylindrical  throat.  The 
"cubic"  nozzle  was  fabricated  in  2",  3"  and  4"  diameters.  Geometry  of  the  cubic  nozzle  is  shown 
fig.29. 

Data  for  the  "cubic"  nozzles  were  measured  repeatedly  over  a  7  year  period  in  both  the  Boeing  Thrust 
Stand  and  the  Boeing  Flight  Simulation  Chamber  (ESC).  The  repeatability  bands  for  the  three  "cubic" 
nozzles  from  both  facilities  are  shown  in  figures  30  &  31.  These  data  were  obtained  when  mass  flow 
accuracy  was  approximately  0.10  to  0.15%,  traceable  to  the  CEESI  300  cubic  foot  Volumetric 
Primary  Standard  (ref.  8).  In  addition,  instmmentation  tolerances  were  left  in  the  data  which  increased 
the  repeatability  bands.  It  should  be  noted  that  the  Velocity  Coefficient  (Cy,  noted  CVN  in  figures  30 
and  31)  for  all  three  nozzle  sizes  is  essentially  the  same,  although  Reynolds  Number  should  increase 
Cy  for  the  4"  nozzle  by  0.1%  above  the  2"  nozzle.  This  increment  can  easily  be  buried  within  the  data 
repeatability.  The  discharge  coefficients  (Cd,  noted  CDN  in  figures  30  and  31)  are  just  as  measured. 
ASME  nozzle  data  repeatability  bands  run  back-to-back  with  the  "cubic"  nozzles,  show  wider  bands 
of  repeatability,  (fig.  32,  33  &  34)  than  the  "cubic"  nozzles.  Both  Cy  and  Cd  show  trends  opposite  to 
normal  Reynolds  Number  effects. 

Due  to  the  questionable  performance  of  the  ASME  nozzles,  both  the  Boeing  Wind  Tunnel  Facilities 
and  the  Boeing  Propulsion  Laboratories  now  use  the  "cubic"  nozzles  to  verify  thrust  facility 
performance. 


42  3"  Cubic  nozzle  test  results  of  different  calibration  tanks 


4.2.1  Model 

The  DNW  3"  cubic  nozzle  is  a  blown  nozzle  which  delivers  thrust  along  the  model  centre  line.  It  is  a 
convergent  nozzle  with  3"  exit  (throat)  diameter.  Its  geometry  is  described  by  a  polynomial  of  the 
third  degree  conformal  to  the  3"  design  used  by  Boeing  (see  fig.29) : 

Y  =  0.046648. X3  -i-  1,5  (inch) 

The  duct  radius  upstream  of  the  nozzle  contraction  is  Rb  =  3.5",  the  nozzle  throat  radius  Rk  =  1.5" 
and  the  nozzle  length  Xmax  =  3.5". 

DNW  owns  two  3"  cubic  nozzles,  identified  as  port  and  starboard  "engine".  A  model  drawing  is 
given  in  figure  35.  The  compressed  air  is  injected  laterally  into  a  settling  chamber  and  passes 
successively  a  throttle  plate  and  "retimet"  mufflers.  The  throttle  plate  orifice  area  is  741  mm2  leading 
to  a  pressure  drop  over  the  plate  of  about  10  bar  at  a  Nozzle  Pressure  Ratio  NPR=1.89  . 

The  nozzle  instrumentation  consists  of : 

-  1  rake  consisting  of  5  total  pressure  probes,  the  pressure  of  the  centre  probe  is  mostly  used  to  set 
datapoint  pressure  ratio  conditions.  The  probes  are  positioned  in  such  a  way  that  an  area  weighting 
procedure  results  by  taking  the  arithmetic  mean  of  the  readings. 
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- 1  rake  consisting  of  3  total  temperature  probes  : 
.  2  Copper/Con  Stan  tan  thermocouples, 

.  1  temperature  resistive  sensor  Pt  100 


4.2.2  Results 

Test  results  on  different  3"  cubic  blown  nozzles  (all  conformal  to  the  3"  design  used  by  Boeing)  are 
provided  by  the  facilities  (calibration  tanks).  The  test  data  are  given  in  the  form  of  nozzle  discharge 
and  velocity  coefficients  (Cd  and  Cy)  as  a  function  of  Nozzle  Pressure  Ratio  (NPR).  The  following 
data  are  available : 

1)  NLR  Engine  Calibration  Facility  :  DNW  starboard  nozzle 

DNW  port  nozzle 

2)  ONERA  S4B  DNW  starboard  nozzle  (two  different  venturis  for  mass  flow 

measurement) 

3)  ARA  Mach  Simulation  Tank  1  ARA  blown  nozzle 

4)  Boeing  Thrust  stand  1  Boeing  blown  nozzle 

In  figures  36  and  37  the  Cd  and  Cy  curves  of  the  DNW  "starboard  nozzle"  are  presented.  Per  facility 
(ONERA  and  NLR)  mean  lines  of  the  test  data  are  shown.  Balance  errors  have  dominating  effects  in 
the  Cy  accuracies.  Therefore  error  bars  have  been  included  in  fig.37,  based  on  0.05%  of  the 
maximum  axial  range  of  the  concerning  balance  systems  (ONERA  S4B  :  3.0N  and  NLR  ECF:  1.5N) 
The  anticipated  errors  become  relatively  large  at  low  NPR  values  (low  thrust)  because  both  calibration 
tanks  are  designed  for  much  higher  loads.  This  explains  the  differences  in  Cy  at  low  NPR. 

The  test  data  of  all  nozzles  (Cd  and  Cy)  are  presented  in  figure  38  and  39  respectively.  Shown  are  the 
mean  lines  of  the  test  data  of  each  facility  and  the  outer  lines  (maximum  and  minimum)  of  all  the 
provided  test  data.  The  same  comments  as  made  above  can  be  made  with  respect  to  the  Cy  curves. 
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4.3.1  Mass  flow 

Boeing's  results  on  the  cubic  nozzle  (fig.30  &  31)  show  that  over  a  long  period  of  time,  it  is  possible 
to  report  mass  flow  measurements  on  choked  nozzles  within  ACd=±0.1%. 

The  absolute  result  should  be  very  close  to  what  is  shown  here  because  of  the  precaution  taken  to 
derive  the  measurements  from  those  of  the  gravimetic  facility  (systematic  error  within  0.07%) 

This  is  to  be  compared  with  comparisons  made  on  the  same  DNW  3"  nozzle  with  three  different 
calibration  tanks.  Although  these  tanks  were  not  calibrated  in  the  same  sophisticated  way  as  Boeing's, 
and  were  using,  either  predicted  values  (ONERA)  or  calibrated  sonic  venturis  (ARA,  NEL)  for  their 
flowmeters,  they  provide  consistent  results.  These  results  differ  by  about  0.015%  and  are  higher  than 
those  of  Boeing  (obtained  on  a  physically  different  nozzle)  by  only  0.08%. 

When  looking  at  aU  data  collected  on  such  nozzles,  the  width  of  the  uncertainty  band  is  about  ±0.2% 
for  choked  cases,  i.e.  about  twice  the  objective. 

4.3.2  Thrust 

Boeing's  data,  (as  well  as  other  data  not  reproduced  here),  fig.30  to  34,  shows  a  good  long  term 
repeatability  within  ACy=±0.2%.  More  interesting  is  to  look  at  what  happens  when  comparisons  are 
made  between  calibration  tanks. 
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The  uncertainty  band  for  choked  conditions  is  about  ±0.3%  overall  and  ±0.17%  if  using  the  same 
physical  nozzle. 

To  summarize,  one  could  say  that  the  state  of  the  art  for  mass  flow  measurement  in  wind  tunnel 
models  is  not  far  from  the  0.1%  objective  for  absolute  accuracy.  Better  understanding  of  nozzle  flow 
conditions  would  help  in  making  this  difficult  goal  more  economically  achievable. 

However,  there  is  a  factor  2  to  3  to  gain  in  absolute  accuracy  on  thrust.  Since  some  wind  tunnels  can 
repeat  tests  with  motorized  models  with  a  drag  repeatability  of  1  count  (0.0001),  it  is  likely  that  this 
can  be  achieved  especially  if  all  precautions  are  taken  as  proposed  in  above  2.3  and  2.4. 


5.  CONCLUSIONS 


-  To  be  consistent  with  the  1  drag  count  (or  0.3%)  accuracy  required  by  wind  tunnel  users,  mass  flow 
and  thrust  of  engine  simulators  must  be  determined  within  ±0.1%. 

-  Survey  of  today's  state  of  the  art  indicates  that  both  goals  could  be  achieved  (for  mass  flow)  or 
approached  (for  thrust)  through  reasonable  improvements  of  existing  methods. 

-  Calibration  by  gravimetry  provides  the  required  accuracy  at  least  at  CEESI  but  limitations  in  mass 
flow  range  make  it  long  and  complex  to  transfer  the  results  to  the  larger  flowmeters  normally  used  in 
wind  tunnels. 

-  Nevertheless,  Boeing  has  developed  a  reliable  methodology  for  such  a  transfer. 

-  Toroidal  throat  nozzles  (short  radius)  must  be  preferred  to  cylindrical  throat  (ASME)  nozzle  for 
sonic  nozzle  flow  measurement. 

-  Several  facilities  seem  today  reasonably  satisfied  with  the  use  of  calculated  mass  flow  for  such 
nozzles,  these  being  partly  validated  through  back  to  back  calibrations  at  different  Reynolds  numbers 
and  interfacility  comparisons. 

-  Nevertheless  progress  should  be  made  in  the  understanding  of  the  boundary  layer  conditions  in 
sonic  nozzles.  Little  has  been  done  on  this  topic  since  1970.  Modem  CFD  is  not  giving  much  more 
reliable  results  yet.  Some  careful,  basic  experiments  would  be  useful.  Such  progress  would  allow  to 
better  design  and  use  sonic  nozzle  flowmeters  and  possibly  to  avoid  gravimetric  calibrations  and  their 
expensive  transfer. 

-  Real  gas  effects  are  properly  accounted  for  with  rules  established  in  the  70's.  Going  further  would 
require  to  check  thermodynamics  tables. 

-  Cubic  nozzles  are  recommended  as  thrust  reference  nozzles  rather  than  ASME  type  nozzles.  Some 
authors  consider  that  subsonic- supersonic  nozzles  giving  a  uniform  flow  in  the  exit  plane  could  be 
even  better  candidates  for  thmst  calibration. 

-  Various  thrust  benches  have  been  presented,  together  with  experimental  data  and  interfacility 
comparisons. 


41 


6.  REFERENCES 


1.  AGARD  CP  n°429  -  Aerodynamic  data  accuracy  and  quality:  requirements  and  capabilities  in  wind  tunnel  testing  - 
Session  III:  drag  accuracy  achievements  (1987) 

2.  J.Leynaert,  Recent  developments  in  wind  tunnel  testing  techniques  at  ONERA  -  ICEFM  1991 
ONERATP  1991-135 

3.  AGARD  AR  n°184  -  Wind  tunnel  flow  quality  and  data  accuracy  requirements,  Nov.  1982. 

4.  Smith  Robert  E.  Jr  and  Matz,  Roy  J."A  Theoretical  Method  of  Determining  Discharge  Coefficients  for  Venturis 
Operating  at  Critical  Flow  Conditions"  ,  Journal  of  Basic  Engineering.  Pages  434  -  446,  Transactions  of  the  ASME, 
Dec.  1962. 

5.  Stratford,  B.  S.,  "The  Calculation  of  the  Discharge  Coefficient  of  IVofiled  Choked  Nozzles  and  the  Optimum  Profile 
for  Absolute  Air  Flow  Measurement",  Journal  of  the  Roval  Aeronautical  Society.  Vol.  68,  pages  237  -  245,  April 
1964. 

6.  Amberg,  B.  T.,  Britton,  C.  L.  and  Seidl,  W.  F. "  Discharge  Coefficient  Correlations  for  Circular  -  Arc  Venturi 
Flowmeters  at  Critical  (Sonic)  Flow",  Journal  of  Fluids  Engineering,  pages  111  -  123,  June  1974. 

7.  Stevens,  Robert  L.,  "Development  and  Calibration  of  the  Boeing  18  kg/s  ( 40  Ibm/sec)  Airflow  Calibration  Transfer 
Standard",  International  Symposium  on  Fluid  Flow  Measurement  Washington  D.C.,  U.S.A.,  Nov.  16  -  19,  1986. 

8.  Arnberg,  B.  T.,  and  Britton,  C.  L.,  "Two  Primary  Methods  of  Proving  Gas  Flow  Meters",  Paper  No.  3  -  8  -  216, 
Symposium  on  Flow,  Pittsburgh,  Pa.,  May  10  -  14,  1971. 

9.  Seidl,  Walter  F.,  "Primary  Calibration  of  the  Boeing  18  Kg/s  (40  Ibm/sec)  Airflow  Calibration  Transfer  Standard”, 
International  Symposium  on  Fluid  Flow  Measurement,  Washington  D.C.,  U.S.A.,  Nov.  16  -  19,  1986. 

10.  Smith,  Robert  E.  Jr.,  and  Matz,  Roy  J.,  "Performance  Characteristics  of  an  8  Inch  Diameter  ASME  Nozzle 
Operating  at  Compressible  and  Incompressible  Conditions",  Presented  at  93rd  Winter  Annual  Meeting,  Nov.  26  -  30, 
1972,  ASME. 

11.  Johnson,  Robert  C.  "Real  -  Gas  Effects  in  Critical  -  Flow  -  Through  Nozzles  and  Tabulated  Thermodynamic 
Properties",  NASA  Technical  Note  NASA  TN  D  -  2565,  Jan.  1965. 

12.  Fromm,  E.  H.,  "The  Boeing  Flight  Simulation  Chamber  for  Static  Calibrations  of  Engine  Simulators",  45th 
Meeting  of  The  Supersonic  Tunnel  Association,  April  1976. 

13. Fluid  Meters,  Their  Theory  and  Application,  Sixth  Edition,  1971,  ASME 

14.  Amberg,  B.T.,  Review  of  Critical  Flowmeters  for  Gas  Flow  Measurements,  ASME,  Journal  of  Basic  Engineering, 
December  1962,  pp  447-460. 

15.  Reimer,  R.M.,  Determination  of  ASME  Nozzle  Flow  Coefficients  by  Thrust  Measurement,  ASME,  Journal  of 
Basic  Engineering,  December  1965,  pp  1058-1062. 

16.  Benedict,  R.P.,  Most  Probable  Discharge  Coefficients  for  ASME  Flow  Nozzles,  ASME,  Journal  of  Basic 
Engineering,  December  1966,  pp  734-744. 

17.  Smith,  R.E.  and  Matz,  R.J.,  Performance  Characteristics  of  an  8-In-Diameter  ASME  Nozzle  Operating  at 
Compressible  and  Incompressible  Conditions,  ASME,  Journal  of  Fluids  Engineering,  December  1973,  pp  542-550. 

18.  Olson,  A.T.,  Nozzle  Discharge  Coefficients  -  Compressible  Flow,  ASME,  Journal  of  Fluid  Engineering,  73- 
WA/FE-7 

19.  Szaniszlo,  A.J.,  Experimental  and  Analytical  Sonic  Nozzle  Discharge 

Coefficients  for  Reynolds  numbers  up  to  8*10^,  NASA  Technical  Memorandum  NASA  TM  X-71595,1974. 


42 


20.  Szaniszlo,  A.J.,  Experimental  and  Analytical  Sonic  Nozzle  Discharge 

Coefficients  for  Reynolds  numbers  up  to  8*  lO^,  ASME,  Journal  of  Engineering  for  power,  74-WA/FM-8. 

21.  Holdhusen,  J.S.,  Some  Recent  Boundary  Layer  Measurements  in  a  5.5  Inch  ASME  Nozzle,  20*^^  Meeting  of  the 
AIAA  Turbine  Engine  Testing  Working  Group,  June  22,  1978,  Minneapolis,  Minnesota,  USA. 

22.  Hilsenrath,  Joseph,  et  al. 

Tables  of  Thermodynamic  and  transport  Properties  of  Air,  Argon,  Carbon  Dioxide,  Carbon 
Monoxide,  Hydrogen,  Nitrogen,  Oxygen,  and  Steam. 

Pergamon  Press,  1960.  New  York. 

23.  B.Masure 

Etude  de  la  precision  des  mesures  de  debits  gazeux  par  la  methode  du  col  sonique. 

Association  Technique  Maritime  et  Aeronautique  session  1968,  ONERA  T.P.  572  (1968). 

24.  Michels,  A.  Wassenaar,  T.  and  Wolkers,  G.  J. 

Thermodynamical  Properties  of  Air  for  Temperatures  between  +75°C  and  -  170°  C  and 
Pressures  up  to  1200  Atmospheres. 

Applied  scientific  Research,  The  Hague,  Section  A.  Vol.  5  (1955)  pp  121  -  136 

25.  John  E.Green 

A  design  of  sonic  nozzle  for  the  precise  measurement  of  mass  flow 
Z.Flugwiss.l9  (1971),  Heft  3 

26.  W.P.Jones  and  B.E.Launder 

The  prediction  of  relaminarization  with  a  two-equation  model  of  turbulence 
InUoumal  of  Heat  and  Mass  Transfer,  Vol.l5,  1972,  pp.301-314 

27.  G.W.Hall 

Application  of  boundary  layer  theory  to  explain  sonic  nozzle  and  venturi  flow  peculiarities 
Proceedings  of  the  institution  of  mechanical  engineers,  Vol. 173,  number  36,  1959 

28.  B.Masure 

Banc  de  dynalpie  au  point  fixe  de  Modane.  Etude  de  I'effet  sur  les  performances  mesurees  d'une  variation  des  conditions 
d'alimentation  de  fair  comprime  a  I’entree  de  la  balance. 

Note  Technique  ONERA  28/414  GAY  (Juin  1972). 

29.  G.Fasso,  M.Pletin,  P.Broussaud 

Installation  pour  I’etude  de  tuyeres  de  reacteurs  a  plusieurs  flux. 

Note  Technique  ONERA  n°  160  (1970). 

30.  B.Masure 

Virial  and  viscosity  effects  on  the  measurement  of  mass  flow  rates  by  means  of  a  sonic  throat  -  Virial  effect  on  nozzle 
thrust. 

Communication  presentee  a  la  28eme  reunion  de  la  "Supersonic  Tunnel  Association"  - 
Denver  -  30  &  3 1  Octobre  1967 

ONERA  T.P.  528  (1967)  et  aussi :  Breve  information  dans:  "Recherche  Aerospatiale  n°120  (sept.oct.1967) 

31.  B.Masure 

Problemes  de  mesure  sur  maquette  de  la  poussee  d'un  arriere-corps  d'avion  supersonique.  Tuyeres  de  reference  ONERA 
T.P.N°  978  (1971)  et: 

Communication  presentee  a  la  reunion  AGARD  sur  les  Entrees  et  Sorties  de  Propulseurs  Aerospatiaux.  Sandefjord 
(Norvege),  13-17  septembre  1971. 

32.  P.Carriere 

Tuyeres  propulsives  -  AGARD  lectures  series  on  supersonic  turbomachinery. 

Varenne  (It^ie),  15  -  20  mai  1967. 

33.  B.Masure 

Private  communication  based  on  ONERA  NT  n°  8/2772  GA,  April  1968  (unpublished). 


43 


List  of  figures 

figure 

Boeing's  transfer  standard  1 

Continuous  curvature  venturi  configuration  2 

Venturi  screening  facility  3 

Repeatability  of  venturi  screening  facility  4 

Throat  diameter  correction  results  5 

CEESI  mass-time  calibration  system  6 

Mass-time  system  calibration  results  for  14  transfer  standard  venturi  meters  7 

NEL  primary  gravimetric  air  flow  standard  -  Line  diagram  8 

General  arrangement  high  pressure  Multiple  Critical  Venturi  (MCV)  9 

MCV  airflow  calibration  Instrumentation  schematic  10 

Zero  thmst  body  for  momentum  tare  evaluation  1 1 

Sketch  of  3"  cubic  nozzle  in  LST;  down-wind  blowing  12  A 

Sketch  of  3"  cubic  nozzle  in  LST;  up-wind  blowing  12  B 

Momentum  tare  evaluation  using  down- wind  and  up- wind  blowing  nozzles  13 

ASME  Long  radius  Nozzle  (LRN)  14 

Reference  condition  for  reference  ASME  nozzle  15 

Numerical  flow  solutions  for  ASME  nozzle  16 

Short  radius  nozzles  17 

Masure's  charts  for  quick  nozzle  flow  and  thrust  prediction  17bis  a,b,c,d 

Numerical  grid  for  ONER  A  short  nozzle  (r=10  mm)  18 

Viscous  flow  solutions  19 

Axial  distribution  of  boundary  layer  quantities  20 

Axial  distribution  of  Launder  and  Jones  parameter  2 1 

Axial  distribution  of  displacement  thickness  22 

Relaminarisation  criteria  23 

Boundary  layer  computations  close  to  the  throat  of  the  ONERA  short  radius  nozzle  24 
Real  gas  effect  on  mass  flow  rate  (Virial  effect)  for  a  sonic  venturi 

Masure’s  results  for  25 

Ratio  of  from  Masure  and  Johnson  26 

Real  gas  effect  on  thrust  (Virial  effect)  for  a  sonic  venturi 

Johnson’s  results  for  Cj.^  27 

Ratio  of  C-i-y  from  Masure  and  Johnson  28 

Boeing's  cubic  nozzle  29 

2"  and  3"  cubic  nozzle  results  -  range  of  uncertainty  30 

4"  cubic  nozzle  results  -  range  of  uncertainty  3 1 

2"  ASME  nozzle  results  -  range  of  uncertainty  32 

3"  ASME  nozzle  results  -  range  of  uncertainty  33 

4"  ASME  nozzle  results  -  range  of  uncertainty  34 

Instmmentation  of  DNW  cubic  nozzle  35 

DNW  3"  starboard  cubic  nozzle  -  Flow  coefficient  Cd  36 

DNW  3"  starboard  cubic  nozzle  -  Velocity  coefficient  Cy  37 

3"  cubic  nozzle  -  Flow  coefficient  Cd  38 

3"  cubic  nozzle  -  Velocity  coefficient  Cv  39 


44 


45 


Figure  1-  Transfer  Standard 


Figure  2-  Continuous  Curvature  Venturi  Configuration 


Figure  4-  Repeatability  of  Venturi  Screening  Facility 


S  {%  OF  MEAN  Cd) 
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NOTE:  SYSTEM  IS  SUBMERGED  IN  A  WATER  8ATH  FROM  SOLENOID  VALVES  TO 
HEAT  EXCHANGER  EXIT 


Figure  €-  Mass-Time  Calibration  System 


Figure  5-  Throat  Diameter  Correction  Results 
( 163  Venturi  Meters) 


r  Cd  =  (-5.03815E-03)(L0Gj(jRd)^  ♦  (6, 39382E-02) (LOGj^Rd)  +  (0.79127) 


Figure 


Mass-Time  System  Calibration  Results  for  14 
Transfer  Standard  Venturi  Meters 
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Fig.  8  -  Line  diagram  of  NEL  primary  gravimetric  air  flow  standard. 
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SECTION  A  -  A 


TOTAL  PERFORATED  AREA 
EQUAL TO  8  TIMES  THROAT  AREA 


High  pressure  Multi  Critical  Venturi  (MCV) 
General  arrangement 


MCV  Code 

Tliroal  Uia. 
(inche) 

1 

0.1870 

2 

0.2645 

4 

0.3740 

8 

0.5289 

16.1 

0.7480 

16.2 

0.7480 

9  i’ 


HIGH  PRESSiiRii; 

MULTIPLE  CRITICAL  VENTURI  -  MCV 


Transfer  Standard 


DEAD  UlEieHr  DEAD  WEIGHT 

TESTER  TESTER 


Figure  10  MCV  airflow  calibration 
Instrumentation  schematic 

MASS  FLOW  TARE  CALIBRATION  ; 

Fcorr  =  ^BAL  “  Kp  (Pz  ~  Fa)  “  Kw 
where : 

Fcorr  =  Balance  measurements  corrected  for  mass  flow  and  pressure  tares 

Fbal  =  Balance  measurements  corrected  for  interactions 

Kp  =  Pressure  tare  constant 

K;v  =  Mass  flow  tare  constant 

Pi  =  Pressure  in  zero  thrust  body 

Pa  =  Ambient  pressure 

Ji  =  Air  temperature  in  zero  thrust  body 

Wa  =  Air  mass  flow 


Pz  , 


Fig.  1 1  -  Zero  thrust  body  for  momentum  tare  evaluation. 
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High  Pressure  Venturi^ 


Fig.  1 2  A  -  Sketch  of  3"  cubic  nozzle  in  LST,  down-wind  blowing.  2  g .  g.  ggj,/g  n^^zje  in  LST,  up-wind  blowing 


Flow 


Down-wind  blowing  nozzle 


Flow 


Up-wind  blowing  nozzle 


Notations  ;  x  :  unit  vector  oriented  down-wind 
T  :  net  thrust  of  the  nozzle 

3  :  horizontal  force  measured  by  the  dynamometer  (and  exerted  on  the  metric  part) 
£  :  momentum  correction  due  to  momentum  effect  at  the  uncoupling  belows. 

Case  :  Tq  .  jr  +  +  '2>f^  -  o 
Cased)  :  Td  .?+e^+3g,  =  o 

If  mass  flow  rates  are  the  same  in  cases  (T)  and  then  |T@|  =  |Tg)|,  strut  deflections 
are  opposite,  so  Tj  ■  a'  =  •  T® .  x  ;  Egj  =  Ejd  =  e. 

Adding  the  two  equalities  gives :  2e  +  3®  +  3®  =  0 
The  momentum  correction  e  can  in  that  way  be  evaluated. 


Fig.  1 3  -  Momentum  tare  evaluation  using  down-wind  and  up-wind  blowing  nozzles. 
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Figure  14  ASME  Long  Radius  Nozzle  (LRN) 


Euler  calculation  (non  viscous) 

•  Air,  Y=  1.4 

•  2D  grid  model 

•  1515  grid  points 

•  2D  axisynunetric  calculation 

Navier-Stokes  calculation  (viscous) 

•  Air,  Y  =  1.4 

•  2D  grid  model 

•  6161  grid  points 

•  2D  axisynunetric  calculation 

•  Turbulent  calculation 

•  Re,  =  8  10‘ 

•  Baldwin  Lomax  turbulence  model 

•  Sutherland’s  formula  for  viscosity 

ll  =  p„(T„+S)/(T  +  S)(T/T„)“ 

S  =  110  K  (air) 


Boundary  conditions 

•  Total  pressure, 

•  Total  temperature. 


Po  =  45  bar 
T„  =  293  K 


UromniSSioSSnSillilli 


Numerical  flow  solutions  for  ASME  nozzle 

a)  Euler  (nonviscous),  reference  nozzle 

b)  Euler  (nonviscous),  with  divergent  exit 

c)  Navier-Stokes  (viscous),  with  divergent  exit 


FIGURE '15  •  Reference  condition  for  reference  ASME  nozzle 


CIRCULAR  ARC  NOZZLE 


SMITH-MATZ  NOZZLE  ONERA  Rj/d;.  NOZZLE 


CONTINUOUS  CURVATURE  ENTRANCE  SHAPE 
(Boeing  Transfer  Standard  Nozzle) 


figure  17  Shon  radius  nozzles 


Figure  1 7  bis  -  a  Mass-flow  rate  calculation 

Sonic  line  curvature  effect  (Masuie  [30]) 


SONIC  LINE 


I  1.106 

l»2  A- 1 


.  Fig.  23  -  Relaminarisation  criteria  (Launder-Jones)  applied  to  ONERA  short  radius  nozzle  ( ^  -  2  ). 

M960312 


Delta  '  (rn)  [displacement  Ihickness] 


Fig.  24  -  Boundary  layer  computations  close  to  the  throat  of  the  ONERA  short  radius  nozzle  ( ^  =  2  ). 
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.-igure  25  Rsai  gas  effecc  or  mass  flow  rate  CYirial  effect)  for  a  sonic  venturi 
Masure's  results  for  Cj^y 
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'igure  26  Real  gas  effect  on  mass  flow  rate  for  a  sonic  venturi 
Rado  of  coefficients  from  Masure  (M)  and  Johnson  (J) 


Figure  27  Real  gas  effect  on  thrust  (Virial  effect)  for  a  sonic  nozzle 
Johnson's  results  for  ^Tv 
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Figure  28  Real  gas  effect  on  thrust  for  a  sonic  nozzle 


Ratio  of  coefficients  from  Masure  (M)  and  Johnson  (J) 


2.0029"  AND  3.0004"  DIA.  ClJlilC  NOZZI.K.S 
1  Iiuus  r  S  I  AND  -  lil.OWlNi;  NOZZ.I.IO  DA  I  A 


Figure  32 


■igure  33 


3.9989"  DIA.  ASME  NO'/ZI-E 


DNW  3"  STARBOARD  CUBIC  NOZZLE 


-5  2.0  2.5  3.0  3.5 

MPR 

Figure  36  DNW  3"  starboard  cubic  nozzle 
Flow  coefficient 
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Fi2ure37  DNW  3"  starboard  cubic  nozzle 
Velocity  coefficient 


Figure  38  3"  cubic  nozzle 
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ON'ERA  static  thrust  benches 
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ONERA  operates  in  the  Modane-Avrieux  Center  three  facilities  devoted  to  measure  nozzles  thrust  and 
mass  flow. 


Two  of  them,  located  in  a  vacuum  tank  named  S4B,  operate  only  with  air  at  ambient  temperature.  The 
third  one,  named  dynalpy  bench  BD2,  can  provide  hot  gas  (up  to  1  lOOK)  and  works  only  at 
atmospheric  exhaust  pressure. 


1.  ONERA  test  benches  located  in  S4B 

The  two  benches  located  in  S4B  are  based  on  the  same  principle  ;  only  components  sizes  are  different 
(balance,  massflowmeters,  plenum  chamber). 

1.1  A  sketch  of  the  facility  which  can  operate  with  a  mass  flow  up  to  3,5  kg/s  is  shown  figure  A) 
(Al.l  &  Ai.2)' 


The  flow  feeding  the  nozzle  is  provided  from  a  high  pressure  storage  (270  bar)  of  very  d^  air  (dew 
Doint  •  below  -60°C).  The  flow,  after  heating  (up  to  80°C)  and  filtering  (12um),  is  controlled  by  a 
regulation  valve  :  the  mass  flow  is  measured  upstream  by  the  way  of  sonic  throats  massflowmeters. 
The  maximum  pressure  upstream  the  massflowmeters  is  9  bar. 

The  airflow  crosses  the  balance  through  an  uncoupling  system  having  a  very  low  stiffness  compared 
to  the  balance  stiffness,  giving  low  and  repeatable  pressure  tares  (<  0,4%  of  loads  capacities),  and 
cancelling  the  momentum  tarel  The  6-components  balance  capacities  are  the  following  : 


X  ;  3300N  Y  :  25000N  Z  :  2()000N 
L  :  3700  Nm  M  :  4000  Nm  N  :  4000  Nm 

The  air  flow  then  feeds  a  metric  plenum  chamber  on  which  the  tested  nozzle  is  connected.  The  6- 
component  balance  directly  measures  the  jet  thrust  vector.  The  pressure  in  the  plenum  chamber  is 
limited  at  a  value  of  1.5  bar  above  the  external  pressure. 


The  facility  is  installed  in  a  sealed  tank  connected  to  vacuum  by  the  way  of  an  adjustable  nozzle, 
controlling  the  exhaust  pressure  in  the  sealed  tank  between  vacuum  and  atmosphenc  pressure. 


Nozzles  with  inlet  at  high  pressure  (blown  nacelle)  (tig.  A.  1.2)  can  be  tested  using  the  primary 
on  the  balance.  The  maximum  pressure  at  the  inlet  of  the  nozzle  is  50  bar  for  a  maximum  mass 

of  1.6  kg/s. 


pipe 

flow 


1  2  A  larger  thrust  bench  is  used  to  test  nozzles  with  a  maximum  mass  flow  of  10  kg/s_.  It  operates  on 
the  same  principle  as  the  previous  one.  The  main  difference  consists  in  the  use  of  a  static  ejector  to  get 
the  pressure  level  at  the  nozzle  exhaust.  The  side  forces  cannot  be  achieved  because  the  l^ge  pressure 
tares  in  these  axis  diston  the  accuracy.  For  this  larger  thrust  bench,  the  capacities  are  the  tollowing. 


Pressure  in  the  plenum  chamber :  1.3  bar  above  external  pressure. 

Thrust  capacity 1000  N  to -1- 5000  N.  ,  ■  .  ■ 

Nozzles  with  inlet  at  high  pressure  (blown  nacelle)  can  be  tested  using  the  primary  pipe  cm  the 
balance.  The  maximum^pressure  at  the  inlet  of  the  nozzle  is  50  bar  for  a  maximum  mass  flow  ot 

4  kg/s. 

The'^configuration  of  the  bench  for  turbofan  simulator  calibration  is  shown  fig.A2. 
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2.  Dynalpy  bench  BD2 

The  dynalpy  bench  BD2  (fig. A. 3)  is  made  up  with  two  frames  (a  fixed  one  and  a  metric  one) 
connected  by  a  dynamometer  measuring  the  axis  thrust  and  by  springy  blades  allowing  to  get  the  side 
forces.  It  can  operate  with  two  airflows  measured  upstream  by  sonic  throats  mass  flow-meters  :  a 
combustion  chamber  on  the  primary  flow  can  heat  a  maximum  mass  flow  of  3,3  kg/s  up  to  1 100  K. 
The  characteristics  of  this  rig  are  the  following  : 

-  thrust  dynamometer  capacity  ;  13500  N  ; 

-  secondary  flow  :  14  kg/s  at  a  maximum  pressure  of  25  bar  ; 

-  primary  flow  :  see  operating  map. 


Figure  A.l  Sketch  of  one  of  the  ONCRA  lest  benches  located  in  S4B 


lijcclor 


Figure  A. 2  Configiiraiion  of  one  of  the  ONFRA  test  benches  located  in  S  ill 
for  turbofan  simulator  calibniiion 


Figure  A. 3  ONERA  dynalpy  bench  BD2 
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APPENDIX  B 
DLR  static  thnist  bench 

The  ETG  is  used  for  calibration  of  Turbine  Powered  Simulators  (TPS)  under  simulated  wind-tunnel 
conditions  without  external  flow.  By  evacuating  the  tank,  a  Mach  number  equivalent  pressure  ratio  is 
simulated  which  is  required  for  determination  of  nozzle  and  thrust  coefficients. 

This  facility  is  also  able  to  measure  nozzle  mass  flow  and  thrust  coefficients. 

A  sketch  of  the  tank  is  shown  fig.Bl;  the  installation  of  a  cubic  nozzle  is  presented  fig.B2. 

Tank  size  :  length  4.2  m  ;  diameter  1.6  m 

Max.  nacelle  diameter :  300  mm 

Tank  pressure  :  .2  x  10^  Pa  up  to  1.5  x  10^  Pa 

Max.  mass  flow  :  up  to  24  k_2/s  at  tank  pressure  1.0  x  10^  Pa 

Vacuum  storage  :  10.000  m^at  .2  x  10^Pa_ 

Drive  air  supply  ;  up  to  4.5  kg/s  at  40  x  10^  Pa 
Drive  air  temperature  :  up  to  415  K. 

Instrumentation 

Six-component  balance  integrated  in  the  front  plate  for  a  maximum  thrust  force  of  6000  N  in  axial 
direction  (accuracy  0.03%). 

9  calibrated  sonic  nozzles  for  high-precision  mass  flow  determination  (accuracy  0.3%). 

Scanivalves,  pressure  transducers  and  PSl  systems  for  a  pressure  range  from  .01  x  10^  up  to  50  x 
105  Pa. 

Statically  and  dynamically  calibrated  thermocouples  and  PT-100  sensors. 

Computer  controlled  4-axls  traversing  mechanism  for  probes  to  measure  the  flow  field  behind  a 
simulator. 

Pressure  and  temperature  rake  for  measurements  upstream  and  downstream  of  a  TPS. 

Data  Acquisition 

Real-time  data  acquisition  and  reduction  with  local  computer  for  numerical  and  graphical  quick  look. 
On-line  data  transfer  to  DLR  network  and  on  request  to  customers. 

Measuring  Techniques 

Measurement  of  thrust,  mass  tlow  rate,  pressure  and  temperature  distribution  behind  fan  and  turbine 
for  determination  of  calibration  coefficients.  Laser  measurements  and  Schlieren  photography  of  the 
jet. 
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Figure  11.1  DLR  Calibialioii  lank 
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Fisure  B.2  Cubic  nozzle  mounted  on  DLR  Calibration  tank 
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APPENDIX  C 

ARA  static  thrust  benches 

1.  The  ARA  Large  Scale  Thrust  Measuring  Rig  LSTMR 

The  LSTMR  facility  (fig. Cl)  provides  a  means  of  assessing  the  performance  of  either  single  or  twin 
stream  nozzle  configurations  both  in  terms  of  nozzle  discharge  and  thrust  characteristics.  A  cylindrical 
pressure  vessel  with  a  central  bulkhead  forms  two  plenum  chambers.  These  two  chambers,  ultimately 
feeding  the  core  and  fan  flows,  tire  connected  through  to  the  ng/model  interface  by  means  of  a  pair  of 
conceiitric  ducts.  Simultaneous  flows  of  up  to  12  kg/s  and  4  kg/s  can  be  provided  for  the  fan  and  core 
respectively.  The  air  supply  to  the  LST.MR  comes  from  two  separate  air  systems.  The  majority  of  the 
flow  is  provided  from  a  .series  of  storage  tanks  initially  pressurised  to  10  bar  with  clean  dry  air.  Air  is 
fed  to  each  chamber  via  a  single  0.250'ni  diameter  air  feed  which  bifurcates  downstream  of  a  master 
control  valve  into  individual  6.25  m  and  0.203  m  lines  for  the  fan  and  core  streams  respectively. 

Other  than  a  pair  of  butterfly  valves  which  govern  the  pressure  drop  into  each  plenum  chamber,  the  10 
bar  supply  is  uncontrolled.  Constant  nozzle  pressure  ratios  (XPRs)  are  achieved  by  augmenting  the 
flow  in  each  stream  with  that  from  a  separate  250  bar  high  quality  air  supply.  The  amount  of 
augmentation  is  controlled  by  a  closed  loop  servo  system  which  responds  to  changes  in  mass  flow. 
The  airflow  in  each  stream  is  ctu'efully  metered  using  critical  venturi  meters.  The  venturi  diameters  in 
the  Fan  and  Core  streams  are  fixed  at  .0089m  and  0.051m  respectively  whereas  those  for  the  250  bar 
augmentation  lines  are  selectable  from  a  minimum  of  0.003m  to  a  maximum  of  0.015m  diameter 
depending  on  the  flow  requirements. 

A  hydrogen  burner  installed  in  the  core  duct  can  be  used  to  provide  heated  air  up  to  600°C  to  the 
nozzle  system  under  test. 

The  plenum  chamber  as.sembly  is  floated  on  a  pair  of  air  bearings  which  constrain  all  loads  with  the 
exception  of  those  in  the  axial  and  torsional  directions.  The  bearings  lift  the  rig  by  0.1  mm  and 
provide  a  near  friction  free  support  enabling  the  thrust  to  be  measured  by  a  single  load  cell  mounted 
on  the  thrust  axis.  The  two  air  bearings  are  mounted  below  the  plenum  chamber  thus  enabling  lateral 
forces  to  be  measured  by  load  cells  mounted  perpendicular  to  the  thrust  axis  at  each  end  of  the  plenum 
chamber  available  with  capacities  ranging  from  2200N  to  8900N. 

2.  The  ARA  Mach  Simulation  Tank  MSTl 

MSTl  was  commissioned  in  1979-80  and  has  subsequently  been  developed  into  a  highly  productive 
and  repeatable  nozzle  system  calibration  facility. 

MSTl,  (fig.C2),  consists  of  a  plenum  chamber  tor  tank)  connected  to  an  external  suction  device, 
providing  tank  pressures  in  the  range  atmospheric  (,M=0)  to  0.35  bar  (.M  =  1.3).  This  variation  (for  a 
fixed  set  of  model  internal  conditions)  is  achieved  by  variation  of  tank  exit  area.  A  binary  set  of  sonic 
venturi  meters  (MCV's)  tu'e  used  to  set  up  an  exit  area  (A.MCV  from  3.23x10-4  m2  to  a  maximum  of 
0.051  m2  in  steps  of  3.23x  lO'^  or  larger).  The  venturis  are  also  used  to  measure  the  tank  exit  mass 
flow  rate. 

The  nozzle  system  under  test  is  mounted  at  the  forwtu-d  end  of  the  tank  supported  on  a  live  (metric) 
frame.  The  forces  acting  on  the  live  frame  and  model  are  measured  by  a  pair  of  six  component  strain 
gauged  balances.  The  balances  are  each  capable  of  measuring  up  to  BOON  axial  force  with  a 
measurement  resolution  of  between  0.22  and  0.44N.  The  large  pressure  area  term  experienced  by  the 
live  frame  due  to  the  differential  pressure  acting  across  it  is  counteracted  by  a  series  of  annular 
pressure  chambers  which  provide  the  reacting  force.  Rolling  diaphragm  seals  are  used  to  maintain 
pressure/mass  flow'  integrit}’  combined  with  a  low  axial  stiftness. 

High  pressure  air  for  TPS.  elector  or  blown  models  is  supplied  from  the  ARA  250  bar  system  capable 
of  providing  air  to  the  model  at  operating  pressures  up  to  69  bar.  The  air  is  transferred  across  the 
balances  by'^ means  of  three  caged  bellows  arrangements  mounted  in  a  plane  perpendicular  to  the 
MSTl  thrust  axis.  A  closed  loop  servo  system  and  high  pressure  sonic  venturi  meter  are  used  to 
control  and  measure  the  mass  flow  to  a  high  degree  of  accuracy. 

Calibrations  using  a  dedicated  system  of  reference  nozzles  tu'e  routinely  performed  to  ensure  the  long 
term  repeatability  of  the  facility. 


^I'OOdiiii 


Figure  C.l  I,:irge  Scale  Tlimsi  Measuring  Rig  (LSTMR)  ai  ARA 
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Figure  C.2  General  arrangement  of  llie  AKA  Macli  Siniulalion  'I'aiik  (MS  I'i ) 
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APPENDIX  D 

NLR  static  thrust  bench 

NLR  Model  Engine  Calibration  Facility  (ECF) 

The  NLR  Model  Engine  Calibration  Facility  (ECF)  is  designed  for  calibration  of  nacelles  equipped 
with  turbo  powered  model  engines  such  as  Turbofan  Propulsion  Simulators.  However,  blown  nacelle 
and  throuah  flow  nacelle  configurations  can  be  calibrated  as  well.  In  the  facility  the  model  engine 
thrust  and  mass  flows  are  detemiined  under  simulated  wind  tunnel  conditions  without  external  flow. 

A  schematical  drawing  of  the  3"  cubic  nozzle  mounted  in  the  facility  is  given  in  fig.D.l.  The  nozzle  is 
mounted  in  the  front  end  of  the  tank.  The  tank  inlet  is  closed  by  the  nozzle  and  the  nozzle  exhausts 
into  the  tank.  Air  suction  from  the  tank  is  established  by  three  centrifugal  compressors  in  series,  with 
a  total  drive  power  of  720  kW. 

A  pipeline  brings  1 10  bar  compres.sed  dry  air  (dew  point  temperature  is  about  218  K  at  the  air  storage 
pressure  of  260  t  r)  to  the  ECF.  In  the  so-called  "skid  ,  the  pressure  can  be  reduced  to  80  bar  and 
below.  The  max;:.', urn  mass  flow  is  6  kg/s.  After  passing  a  heat  exchanger  (heating-up  the  drive  air  to 
a  temperature  of  293-343  K),  the  compressed  air  has  to  be  supplied  to  the  model  through  an  air  duct 
which  has  to  cross  the  balance  with  only  minor  interaction.  Three  flexible  air  couplings  are  mounted, 
each  with  two  degrees  of  freedom,  so  that  balance  forces  exerted  by  the  through-flowing  high 
pressure  air  in  all  six  possible  degrees  of  freedom  are  minimized.  The  high  pressure  air  mass  flow  is 
measured  in  the  supply  line  by  a  sonic  venturi. 

The  mass  flow  can  also  be  measured  at  the  downstream  end  of  the  tank  via  a  set  of  9  sonic  venturis. 
When  they  are  choked  their  mass  flow  can  be  detemiined  ver\'  accurately.  At  steady  state  conditions 
the  hitih  pressure  mass  flow  should  be  equal  to  the  tank  venturis  mass  flow.  Because  steady  state  is 
approached  asymptotically,  the  tank  venturis  mass  flow  is  used  for  check  only.  During  a  test 
sufficient  tank  venturis  are  opened  to  run  them  unchoked  during  most  of  the  test  datapoints.  The  tank 
pressure  can  then  be  controlled  independently  by  the  suction  compressors  rpm.  Only  at  the  highest 
mass  flows  sonic  conditions  are  attained  tor  a  mass  flow  check. 

For  measurement  of  the  model  engine  thrust  the  tiicihty  is  equipped  with  an  e.xtemal  balance  system. 
The  balance  consists  of  bars  which  connect  the  metric  part  to  earth.  Three  bars  are  connected  through 
high  accuracy  load  cells  so  that  three  components  are  measured  (side  force,  axial  force  and  yaw 
moment). 

A  central  bellows  forms  a  flexible  air-tight  connection  between  the  metric  part  of  the  balance  and  the 
non-metric  tank.  In  this  central  bellows  the  engine  model  is  mounted  60  mm  below  the  centre  line  (see 
fig.Dl).  To  minimize  the  parasitic  force  due  to  the  pressure  differences  over  the  central  bellows  two 
compensating  bellows  are  installed  at  the  upper  and  lower  side  of  the  central  bellows. 

The  main  characteristics  are  the  following  : 

-  internal  diameter  of  the  nacelle  housing:  ().5m 

-  maximum  engine  weight  (incl.support  ):  15()kg  at  ().8m  from  balance  center 

-  maximum  axial  thrust  force:  :r3,()()()N 

-  maximum  side  force:  3:900N 

-  fan  inlet  air:  atmospheric 

-  drive  air:  maximum  6  kg/'s  at  80  bar 

-  drive  air  control:  incremental  steps  of  0.0005  kg/s 

-  drive  air  temperature:  243  K  to  343  K 

-  maximum  tank  suction  capacity:  0.58  bar  (M=0.92) 

-  tank  size:  L=6M  D=3M 

-  data  processing  and  data  plotting:  on  line 

-  fan  exhaust  mass  flow:  see  map 


69 


APPENDIX  E 

Boeing  Thnist  Stands  or  Benches 

Four  static  thnist  stands  are  provided  at  Boeing's  Seattle,  Washington  plant  for  thrust 
and  airflow  measurements.  Two  thrust  stands  are  at  the  Wind  Tunnel  Facilities  and  two  thrust 
stands  are  at  the  Propulsion  Laboratories.  Each  thrust  stand  has  a  unique  capability  and 
purpose. 

Two  continuous  air  supplies  are  shared  by  the  Wind  Tunnel  Facilities  and  the 
Propulsion  Laboratories.  A  high  pressure  system  supplies  9  kg/s  at  69  bar.  A  low  pressure 
system  supplies  20.4  kg/s  at  21^  bar.  The  air  is  filtered  and  dried  to  a  dew  point  below  -  46°C. 
Both  systems  have  storage  so  the  above  continuous  flow  rates  can  be  exceeded  for  limited  time 
periods. 

2-Balance  Thrust  Stand 

Tnc  fust  of  the  Wind  Tunnel  thrust  benches  is  the  2-Balance  Thrust  Stand  shown  in 
figure  E.l.  Loads  are  measured  by  two  6-component  balances  mounted  1.829  rneters  apart.  The 
main  load  carrying  member  is  located  0.457  meters  below  the  balance  centerlines,  so  that  the 
thrust  vector  can  be  located  at  or  near  the  plane  of  the  balance  axial  force  centerlines  as  shown. 
A  complete  wind  tunnel  model  with  engine  nozzles  on  both  wings  can  be  calibrated  with  the 
thrust  vectors  on  or  close  to  the  plane  of  axial  force.  The  nozzles  can  be  calibrated  separately  or 
together.  The  type  of  air  flexure  shown  here  has  no  momentum  tare  and  only  a  small  repeatable 
pressure  tare.  The  balances  are  thermally  isolated  from  temperature  gradients  by  insulators  at 
both  ends  of  each  balance.  Load  induced  and  thermally  induced  dimensional  changes  of  the 
stmctures  between  the  balances,  produce  equal  and  opposite  loads  in  the  balances  which  cancel. 
The  twelve  components  measured  by  the  two  balances  are  resolved  into  6-component  data 
midway  between  the  two  balances.  Thrust  loads  up  to  2900  N  can  be  measured,  accurately 
defining  the  magnitude,  angularity,  and  position  of  the  thrust  vector.  Alternate  balances  can 
increase  thrust  capacity  up  to  5300  N.  Air  mass  flow  can  be  measured  upstream  of  the  thmst 
stand  with  a  Multiple  Critical  Venturi  (MCV)  shown  in  Figure  9,  or  with  sonic  venturis 
mounted  inside  the  wind  tunnel  model.  Total  model  weight  can  be  up  to  550  kg. 

Flight  Simulation  Chamber  (FSCi 

The  second  Wind  Tunnel  thrust  facility  is  the  Right  Simulation  Chamber_(FSC)  which 
was  designed  primarily  for  calibrating  Turbo-Powered-Simulators  (TPS).  This  facility  is  shown 
in  Figure  E.2  and  described  in  more  detail  in  Reference  12.  In  addition  to  TPS  calibrations,  it  is 
used  to  calibrate  blowing  nacelles  and  define  the  internal  drag  of  flow-through  nacelles.  Like  the 
2-Balance  Thrust  Stand,  loads  are  measured  by  two  6-component  balances,  and  the  same  type 
of  air  bridge  is  used  for  high  pressure  air.  Thrust  can  be  measured  up  to  1550  N. 

High  pressure  air  flow  can  be  measured  upstream  of  the  air  bridge  by  a  sonic  venturi 
or  high  pressure  MCV.  The  total  flow  through  the  chamber  is  measured  by  the  Low  Pressure 
MCV  at  the  downstream  end  of  the  chamber.  Air  mass  flow  ranges  for  High  Pressure  and  Low 
Pressure  MCV's  are  shown  in  Figure  E.3.  Two  air  ejectors  provide  suction  to  maintain  sonic 
flow  in  the  Low  Pressure  MCV.  Model  weights  up  to  180  kg  can  be  accommodated. 

A  7  year  history  of  standard  "cubic"  thrust  nozzle  data  from  the  2-Balance  Thmst 
Stand  and  the  FSC  are  shown  in  figures  30  and  31. 
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Large  Dual  Flow  Rig  (LDFR) 

The  Propulsion  Laboratory  Large  Dual  Row  Rig  (LDFR)  is  designed  for  static  nozzle 
thrust  and  thrust  reverser  testing.  Loads  are  measured  by  a  3-component  b^ance  in  a  horizontal 
plane,  with  thrust  capacity  from  -  2225  to  8900  N.  Two  major  air  bridges  bring  primary  and 
secondary  airflow  from  the  low  pressure  system  on  to  the  balance.  The  primary  airflow  can  be 
heated  on  balance  by  a  Propane  burner  system  to  temperatures  up  to  922°K  with  flow  rates  up  to 
9  kg/s  and  pressure  up  to  9.3  bar.  The  secondary  flow  system  is  ambient  temperature  at  flow 
rates  up  to  15.9  kg/s  and  pressure  up  to  9.3  bar.  The  primary  and  secondary  flow  systems  are 
independently  controlled,  and  flow  rates  are  measured  by  large  MCV's  upstream  of  the  balance. 
A  tertiary  flow  system  has  been  added  to  provide  air  flow  up  to  2.7  kg/s  and  pressure  up  to  48 
bar  for  complete  simulation  of  some  nozzle  systems.  A  third  pressure  control,  MCV,  and  air 
bridge  is  provided  for  this  system.  Model  weight  up  to  725  kg  can  be  accommodated. 


High  Pressure  Ratio  Rig  ('HPR') 

The  second  Propulsion  Laboratory  thrust  bench  is  the  High  Pressure  Ratio  Rig  (HPR) 
which  can  provide  a  single  flow  at  nozzle  pressure  ratios  from  2  to  40.  Loads  are  measured  by  a 
6-component  balance  with  thrust  capacity  up  to  ±  4450  N.  The  single  flow  system  can  provide 
air  mass  flow  up  to  11.3  kg/s  at  pressures  up  to  48  bar.  Air  mass  flow  is  measured  by  a  high 
pressure  MCV.  Model  weight  up  to  725  kg  can  be  accommodated. 

Each  of  the  Propulsion  Laboratory  thrust  rigs  is  housed  in  a  6  by  6  meter  cross  section 
by  19  meter  long  test  cell  vented  at  the  far  end  to  maintain  atmospheric  pressure  in  the  cell,  and 
isolate  noise  from  the  surrounding  areas. 


2  -  nAlANCF.  TimilST  STAND 


Figure  E. 1 
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APPENDIX  F 


ONERA  SHORT  RADIUS  NOZZLE  (Rj,/dc  =  2) 


1-  Notations  and  conventions 


1-1  Axis  system 

The  selected  axis  system  is  Oxy;  the 
point  O  is  the  geometrical  center  of  the 
throat;  Ox  axis  is  the  throat  axis,  positive 
from  upstream  to  downstream;  Oy  is  ortho¬ 
gonal  to  Ox,  positive  to  the  top. 

1-2  Throat  definition 

The  throat  shape  is  broken  up  in  three 
parts  :  from  downstream  to  upstream: 

-  a  circular  shape;  the  center  is  the 
point  P,  and  the  radius  R^,  is  equal  to  two 

times  the  throat  diameter  dc.  R^.  =  2d(, 

This  part  stops  at  the  point  M,  at  an  angle 
of  60°.  The  coordinates  of  P  and  M  are  the 
following: 


r 


p 


xp  =  0 

yp  =  --de 
2 


M 


xm  =  -f3  dc 


-  a  linear  part,  from  M  to  B,  which  is  a  part 
of  the  tangent  in  M  to  the  circle  (P,  Rc);  the 
tangent  equation  is: 


y  =  13  X  +  2-dc 
7 


-  a  second  circular  curvature;  the  cen-  » 
ter  is  the  point  T;  the  radius  is  r.  C  and  B  ^ 


are  on  the  circle;  BM  is  the  tangent  in  B  to  the  circle  (T,  r). 
the  coordinates  of  C  are  given  values: 


C 


xc  =  -  58,5  mm 
yc  =  -  51,2  mm 


2-  Coordinates  of  points  T  and  B  and  value  of  radius  r 

2-1  Ideal  case 

The  ideal  case  is  to  get  an  horizontal  tangent  in  C.  In  this  case,  the  coordinates  of  the  different 
points  would  be: 


T 


XT  =  Xc  =  -58,5  mm 
yT  =  -VTxc  +  2yc-l,5dc 


B 


xb  =  -  -  xc  +  ^yc  -  ^^dc 
2  2  4 


ye 


=  -  ^xc  +  -yc  -  -d< 
2  2  4 


r  =yc  -  f^xc  - 1,5  dc 


2-2  Real  case 

In  fact,  for  machining  problems,  it  is  difficult  to  get  an  horizontal  tangent  in  C.  So,  we  give  for  ra¬ 
dius  r,  a  value  slightly  greater  than  the  one  defined  in  2-1.  Once  the  value  of  r  is  selected,  the  coordi¬ 
nates  of  T  et  B  can  be  calculated. 
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